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Abstract 


The  subject  of  this  thesis  is  the  system  consisting 
of  the  engagements  between  attacking  intercontinental  ballis¬ 
tic  missile  (ICBM)  reentry  vehicles  (RV)  and  a  defending  bal¬ 
listic  missile  defense  system.  The  thesis  present::  a  brief 
overview  of  the  actual  system  before  proceeding  with  the 
development  of  a  computer  simulation  model  designed  to  aid 
analysis  of  the  system.  The  primary  language  of  the  main 
program  is  SLAM;  the  supporting  programs  use  FORTRAN  V.  The 
RV/Antiballistic  missile  (ABM)  system  is  modeled  as  a  net¬ 
work  through  which  the  RVs  flow.  The  capabilities  of  the  RV 
with  regard  to  yield  and  accuracy  can  be  set  to  the  user’s 
desires  as  can  the  vulnerabilities  of  the  RV  to  the  two  pri¬ 
mary  ki'l  mechanisms  considered.  X-ray  and  neutron  radiation. 
The  ABM  system  consists  of  a  two-tiered,  layered  terminal 
defense  system  with  high  altitude  (70,000  -  150,0.d0  feet) 
and  low  altitude  (10,000  -  70,000  feet)  interceptors.  Inter¬ 
ceptor  yield  and  accuracy  are  changeable.  Either  layer  can 
be  activated  or  deactivated,  singly  or  together,  to  permit 
flexibility  in  the  comparison  of  actual  or  hypothetical  sys¬ 
tems.  An  example  simulation  is  accomplished  to  demonstrate 
model  operation  and  permit  systemic  analysis.  The  thesis 
also  contains  a  User’s  Manual  for  those  interested  in  using 


the  model. 


A  SIMULATION  MODEL  FOR  ANALYZING 
REENTRY  VEHICLE/ ANTIBALLISTIC  MISSILE  ENGAGEMENTS 


I  Introduction 


General  Situation  and  Importance 

In  their  1S80  report  on  military  posture  to  Congress, 
the  Organization  of  the  Joint  Chiefs  of  Staff  collectively 
stated: 

The  intercontinental  ballistic  missile  (ICBM)  force 
provides  flexibility  for  employment  throughout  the  spec¬ 
trum  of  nuclear  conflict— from  limited  options  to  full- 
scale  retaliation.  The  combination  of  responsiveness, 
accuracy,  reliability,  rapid  retargeting  capability, 
and  assured  penetration  makes  the  land-based  missile 
well  suited  for  attack  options  where  promptness,  mini¬ 
mum  collateral  damage,  and  high  confidence  of  success 
are  paramount.  The  flexibility  of  employment  which  the 
ICBM  force  offers  the  National  Command  Authority  is  cf 
great  importance  in  nuclear  planning  (Ref  10:41). 

In  addition  to  the  capabilities  listed  in  the  state¬ 
ments  above,  the  ICBM  force  has  unique  characteristics  which 
contribute  to  the  vital  diversity  of  the  strategic  Triad  of 
bombers,  ICBMs,  and  sea-launched  ballistic  missiles  (SLBM). 
These  unique  characteristics  include  the  Soviet  perception 
of  ICBM  importance,  high  alert  rates,  low  operating  costs, 
a  time-critical  hard  target  capability,  and  defense  suppres¬ 
sion  (Ref  24:15). 

In  considering  the  importance  of  the  ICBM  contribu¬ 
tions  to  the  Triad  and  strategic  stability,  the  perceptions 
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of  the  Soviets  must  be  considered.  Traditionally,  they  have 
concentrated  the  great  majority  of  their  strategic  firepower 
in  their  ICBM  force.  Their  massive  research  and  development 
efforts  and  deployment  plans,  including  four  new  ICBMs,  give 
dramatic  evidence  of  the  importance  they  continue  to  attach 
to  the  ICBM  force.  The  Strategic  Arms  Limitation  negotia¬ 
tions  of  the  1970s  granted  the  Soviets’  ICBM  numerical  and 
size  superiority.  They  clearly  feel  that  the  ICBM  is  vitally 
important  in  maintaining  perceived  essential  equivalence 
(Ref  24:15). 

The  high  alert  rates  of  the  ICBM  force  are  another 
clear  Indicator  of  the  importance  of  the  ICBM.  Normally, 
almost  100  percent  of  the  U.S.  ICBM  force  is  in  instant  readi¬ 
ness  for  launch,  compared  with  about  one-fourth  of  the  bomber 
force,  and  somewhat  less  than  50  percent  of  the  SLBMs .  The 
fact  that  this  high  alert  rate  is  maintained  with  only  about 
one-eighth  of  the  total  operating  budget  of  the  entire  Triad 
also  buttresses  the  case  for  the  significance  of  the  ICBM 
(Ref  24:15). 

Because  of  their  accuracy  and  instant  readiness, 

ICBMs  have  a  time-critical  hard  target  capability,  as  well. 
Combined  with  their  short  time  of  flight,  the  ICBMs  can  both 
limit  the  destruction  of  this  country  in  a  full-scale  nuclear 
exchange,  by  destroying  enemy  forces  before  they  can  be 
launched,  and  act  as  a  force  multiplier  for  other  parts  of 
the  Triad.  The  ICBM  can  be  used  in  a  defense  suppression 
role,  thereby  enhancing  the  effectiveness  of  the  bomber  force. 


In  summary,  the  ICBM  offers  vital  contributions  to 
the  strategic  Triad  of  forces  and  the  maintenance  of  politi¬ 
cal  stability.  But  the  survivability  of  the  ICBM  has 
recently  been  questioned.  Many  experts  feel  that  the  fixed 
locations  of  the  ICBMs  make  them  vulnerable  to  attack.  This 
lack  of  prelaunch  survivability  has  a  severe  impact  on  the 
credibility  of  the  ICBM  deterrent.  Correcting  this  defi¬ 
ciency  continues  to  be  the  highest  priority  strategic  initia¬ 
tive  (Ref  24:3)-  The  billions  of  dollars  that  will  be  spent 
to  make  the  MX  survivable  will  be  wasted,  however,  unless 
the  survivability  of  the  ICBM  is  ensured  during  penetration 
of  enemy  defenses.  There  are  two  sides  to  ICBM  credibility: 
the  ICBM  must  be  capable  of  withstanding  an  enemy  attack,  and 
it  must  be  capable  of  striking  enemy  targets  with  a  high  prob¬ 
ability  of  success.  The  U.3.  government  is  making  massive 
expenditures  to  ensure  the  first  aspect  but  the  second  is  in 
doubt . 

An  effective  ballistic  missile  defense  (BMD)  system 
can  severely  degrade  or  neutralize  the  deterrent,  capability 
of  the  ICBM  force.  Opponents  of  the  Safeguard/Sprint  BMD 
systems  of  the  1960s  characterized  the  defense  as  "trying  to 
hit  a  bullet  with  a  bullet."  Technological  advances  have 
invalidated  this  metaphor.  The  Soviets  clearly  believe  that 
the  technology  is  at  hand  to  deploy  a  credible  and  effective 
ballistic  missile  system  (Ref  6). 

Test  and  installation  of  improved  BMD  systems  are  now 
prohibited  by  the  Anti-Ballistic  Missile  (ABM)  Treaty  of  1972 
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As  a  result,  the  level  of  U.S.  concern  about  BMD  has  fluctu¬ 
ated.  An  opportunity  exists  for  ABM  treaty  reconsiderations 
in  1982.  At  that  time,  the  United  States  and/or  the  Soviet 
Union  can  propose  modifications  to  or  continue  with  the 
treaty. 

The  Soviets  could  also  decide  to  install  an  improved 
BMD  system  unilaterally.  That  is,  they  could  simply  refuse 
to  engage  in  treaty  negotiations,  allow  the  current  treaty 
to  lapse,  and  begin  deployment  of  an  improved  ABM  system. 

Since  deterrence  relies  on  the  perceptions  of  an  enemy,  it 
is  only  necessary  for  the  enemy  to  believe  that  their  defen¬ 
sive  network  is  effective  for  deterrence  to  fail.  As  a 
result,  the  deployment  of  such  a  system  could  significantly 
affect  the  credibility  of  the  United  States'  ICBM  deterrent 
force  and  could  result  in  the  ultimate  in  political  instabil¬ 
ity,  the  failure  of  deterrence  and  the  occurrence  of  nuclear 
war. 

The  analysis  of  the  impact  that  such  a  deployment 
would  have  is  an  extremely  complex  problem.  Analytical  solu¬ 
tions  do  not  exist  and  current  computer  simulation  models 
are  too  large  and  unwieldy,  not  comprehensive  enough,  or 
biased  by  the  designer's  frame  of  reference. 

Literature  Review 

As  illustrated  by  Figure  1,  a  wealth  of  analytical 
approaches  existed  as  early  as  1968  to  aid  the  planner  in 
solving  the  general  problem  of  allocating  missiles  to  targets. 
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The  so-called  missile  allocation  problem  has  four  discrete 
parts  according  to  Matlin:  the  weapon  complex,  the  target 
complex,  the  engagement  model,  and  the  damage  model  (Ref  12: 
338). 

The  weapon  complex  portion  of  the  problem  can  be  de¬ 
scribed  by  three  characteristics:  scope,  reach,  and  commit¬ 
ment.  The  scope  of  the  weapon  complex  refers  to  the  number 
of  different  RV  types  considered  and  whether  or  not  penetra¬ 
tion  aids  are  considered.  Weapon  reach  refers  to  which  weapon 
can  reach  which  targets  and  with  what  degradations  in  accu¬ 
racy  and  payload.  Weapon  commitment  policy  indicates  the 
number  of  waves  launched,  quality  of  damage  assessment,  and 
weapon  availability  uncertainties.  The  target  complex  is 
also  characterized  in  three  categories.  The  scope  of  the 
target  refers  to  dependence  or  independence  of  point,  area, 
and  collateral  targets.  Target  value  defines  an  appropriate 
measure  of  effectiveness  of  a  given  missile  allocation.  The 
target  defense  type  may  be  terminal,  area,  or  preferential 
(Ref  12:337). 

Figure  1  represents  36  different  analytic  approaches 
that  existed  in  1968.  However,  none  of  the  models  represented 
were  comprehensive  in  all  areas.  Generally,  if  the  weapon 
complex  portion  of  the  model  was  sophisticated,  the  target 
complex  or  the  engagement  model  was  simplified  (Ref  12:337). 

In  addition,  these  analytic  techniques  are  aimed  toward  ana¬ 
lyzing  a  different  problem  than  determining  the  impact  of 
ABM  deployment.  The  complexities  and  stochasticities  of  the 


problem  are  such  that  no  analytical  technique  exists  to 
approach  the  problem. 


In  an  effort  to  remedy  the  shortcomings  of  the  anc- 
lytical  approaches,  some  analysts  turned  to  computer  simula¬ 
tion.  The  U.S.  Air  Force’s  Multiple  Engagement  Model  and  the 
McDonnell  Douglas  Astronautics  Company’s  ENGAGE  model  are 
examples  of  some  current  efforts  in  the  area  (Ref  22 j  26). 
Current  simulation  models  do  not,  however,  adequately  analyze 
the  impact  of  ABM  deployment.  Some  are  not  comprehensive  in 
all  areas  of  the  problem.  Those  that  embody  most  of  the  com¬ 
plexities  are  too  large  and  unwieldy.  Some  models  are  also 
biased  by  the  designer’s  frame  of  reference  or  parochialism. 

This  is  not  intended  to  degrade  the  efforts  of  pre¬ 
vious  modelers.  On  the  contrary,  their  models  were  and  are 
worthwhile  because,  in  the  absence  of  analytic  techniques, 
computer  simulation  offers  the  best  approach  to  the  problem. 
In  the  first  place,  since  any  postulated  BMD  system  must  be 
hypothetical,  analysis  by  experimentation  with  the  actual 
system  is  not  possible  nor  practical.  Secondly,  the  process 
of  development  of  a  simulation  model  leads,  by  itself,  to  a 
fuller  understanding  of  the  actual  system  and,  thus,  better 
analysis.  Finally,  the  complexities  and  uncertainties  of 
the  actual  system  have  prevented  the  development  of  analyti¬ 
cal  solutions.  As  an  example  of  the  complexity  of  the  issues 
involved,  consider  the  deployment  of  the  Advanced  Maneuvering 
Reentry  Vehicle  (AMARV). 

The  AMARV,  unlike  current  reentry  vehicles  (RV)  that 
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essentially  freefall  to  the  target  area,  can  maneuver  after 
separation  from  the  ICBM  that  delivered  it.  This  maneuvering 
capability  theoretically  increases  the  possibility  of  success¬ 
ful  penetration  of  an  ABM  defense  network.  In  order  to  ma¬ 
neuver,  guidance  equipment  is  required  inside  the  AMARV  which 
displaces  some  of  the  nuclear  material  that  could  be  carried 
on  board  a  nonmaneuvering  RV,  resulting  in  smaller  weapon 
yields  (Ref  1:43-45).  The  complex  trade-off  question  is  then: 
Should  the  AMARV  be  deployed  in  light  of  its  increased  prob¬ 
ability  of  successfully  penetrating  ballistic  missile  defense 
systems  or  will  its  smaller  yield,  when  compared  to  current 
RVs,  result  in  less  target  destruction  regardless  of  its  pene¬ 
trability?  Current  techniques,  both  analytical  and  using 
computer  simulation,  are  simply  incapable  of  supplying  an 
answer  to  questions  such  as  this.  The  model  contained  in 
this  thesis  is  specifically  designed  to  rapidly  provide  infor¬ 
mation  to  help  answer  such  difficult  questions. 

Problem  Statement 

The  system  under  consideration  consists  of  three 
major  components:  the  attacking  ICBMs  with  their  RVs  and 
penetration  aids  such  as  decoys  and  chaff,  the  defending  ABM 
subsystem,  and  the  target  area.  In  order  to  analyze  the 
effects  of  improving  the  defenses,  the  complexities  of  the 
entire  system  must  be  embodied  in  a  model.  For  example, 
parameters  affecting  the  defense,  such  as  weapon  yields  and 
accuracies,  intercept  altitudes,  and  command  and  control 
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system  reliabilities,  must  be  included.  The  performance 
requirements  of  the  model  can  also  be  developed  in  light  of 
the  questions  the  model  can  be  used  to  answer.  Along  with 
the  previous  question  concerning  AMARV  deployment,  alterna¬ 
tive  RV  designs  can  be  evaluated,  competitive  defense  systems 
and  strategies  can  be  compared,  or  new  penetration  aids 
tested.  In  light  of  the  flexibility  required  for  a  model 
that  must  provide  information  to  help  answer  questions  such 
as  these,  and  the  extreme  complexity  of  the  actual  system, 
the  following  problem  statement  is  appropriate: 

Develop  a  computer  model  of  the  RV/ABM  engagement 
system  that  has  the  flexibility  to  vary  the  param¬ 
eters  affecting  the  capabilities  of  the  three 
major  components  of  the  system:  the  attacking 
reentry  vehicles  and  penetration  aids,  the  defend¬ 
ing  ballistic  missile  defense  network,  and  the  tar¬ 
get  area.  Combined  with  this  flexibility  must  be 
the  sophistication  necessary  to  analyze  the  complex 
interactions  between  the  participants  in  all  parts 
of  the  system. 

This  problem  statement  leads  to  the  objectives  of  this  thesis 
as  outlined  in  the  next  section. 


Objectives 

Claims  for  RV  and  ABM  capabilities  are  advanced  by 
all  parties  to  the  debate.  The  first  objective  of  this  the¬ 
sis  is  to  provide  a  model  that  allows  the  comparison  of  dif¬ 
ferent  RVs,  ABMs,  and  defense  and  offense  strategies.  The 
key  word  for  this  model  is  flexibility.  It  must  be  able  to 
handle  variation  in  parameters,  capabilities,  ana  strategies. 
The  second  major  objective  of  this  thesis  is  a  User’s  Manual 
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that  will  provide  an  indepth  explanation  of  the  model  and 
instructions  on  operation  for  anyone  who  desires  to  use  the 
model.  The  final  major  objective  is  an  illustrative  exercise 
of  the  model  that  demonstrates  both  the  capabilities  of  the 
model  and  the  interactions  of  the  actual  system.  This  exam¬ 
ple  also  permits  the  development  of  some  conclusions  concern¬ 
ing  the  system. 

The  Engagement  Regime 

The  system  that  is  the  subject  of  this  thesis  involves 
the  encounters  between  attacking  ICBMs,  and  the  weapons  car¬ 
ried,  and  defending  ABMs.  These  engagements  take  place  in 
the  atmosphere  above  the  target  area.  The  ICBMs  deploy  indi¬ 
vidual  reentry  vehicles  (RV)  which  enter  the  atmosphere  along 
with  the  chaff  and  decoys  deployed  with  them  to  deceive  the 
defensive  radars.  The  ABM  radars  can  begin  to  discriminate 
the  threat  objects  as  to  actual  type  after  the  atmospheric 
drag  has  stripped  the  chaff  screen  away.  Final  identifica¬ 
tion  of  threat  type  (RV  or  decoy)  takes  place  nearer  the  tar¬ 
get  (Ref  13:12-13).  Figure  2  shows  an  example  altitude  of 
150,000  feet  for  this  final  identification  point.  This  point 
also  marks  the  upper  boundary  of  the  endo-interceptor  region 
or  that  region  where  attacking  RVs  are  engaged  only  by  endo- 
lnterceptors.  The  lower  boundary  of  this  region  is  assumed 
to  be  the  minimum  altitude  at  which  the  endo-interceptor  can 
successfully  engage  the  attacker.  Figure  2  illustrates  this 
boundary  at  70,000  feet  (Ref  6). 


Current  ABM  interceptors  employ  the  effects  of  nuclear 
weapons  to  destroy  the  attacking  RVs.  These  effects  include 
X-rays,  neutron  radiation,  and  blast.  Only  the  first  two 
effects  are  of  primary  concern  because  they  have  the  far¬ 
thest  lethal  radii  (Ref  4). 

In  a  vacuum.  X-ray  radiation  is  lethal  at  a  greater 
distance  than  neutron  radiation,  all  other  factors  being  equal. 
However,  X-rays  are  more  quickly  attenuated  by  the  atmosphere 
than  neutrons.  As  a  result,  the  primary  kill  mechanism 
throughout  the  majority  of  the  endo-intercept  (El)  region  is 
X-ray  radiation  because  of  the  low  density  of  the  atomosphere. 
The  point  intercept  (PI)  region,  on  the  other  hand,  has 
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denser  air  and  the  dominant  kill  mechanism  is  neutron  radia¬ 
tion  or  fluence. 

The  top  of  the  PI  region  is  marked  by  the  maximum 
operational  altitude  of  the  TT  interceptor.  Intercepts  can 
take  place  anywhere  in  the  region.  The  lowest  point  of  the 
PI  region  is  the  minimum  engagement  altitude  of  the  inter¬ 
ceptor  and  results  in  the  Dead  Zone  depicted  in  Figure  2, 
a  region  of  no  engagements. 

Methodology 

The  research  methodology  employed  in  this  thesis  falls 
within  the  framework  of  the  system  science  paradigm  as  out¬ 
lined  in  Management  Systems  -  Conceptual  Considerations  (Ref 
20:295-305).  The  first  step  is  the  conceptualization  of  the 
system.  In  the  conceptualization  of  the  system,  the  rela¬ 
tionships  between  the  system  components  and  how  they  inter¬ 
act  are  explored.  A  critical  step  in  this  phase  is  the  devel¬ 
opment  of  a  structural  model  of  the  system. 

The  second  phase  of  the  system  science  paradigm  is 
the  analysis  and  measurement  phase.  In  this  portion  of  the 
methodology,  such  topics  as  the  selection  of  a  computer  lan¬ 
guage,  the  specifications  of  the  model,  and  the  standardiza¬ 
tion  of  the  model  results  are  dealt  with.  This  definition 
of  the  desired  output  of  the  model  leads  to  the  third  phase 
of  the  paradigm,  computerization. 

The  computerization  phase  of  the  methodology  deals 
with  the  translation  of  the  concepts  of  the  actual  system 
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(the  real  world)  and  the  structural  model  into  a  computer 
model  that  will  satisfy  the  first  objective  of  this  thesis. 

Overview 

After  providing  a  brief  explanation  and  description 
of  the  actual  system  in  Chapter  II,  the  remainder  cf  this 
thesis  deals  with  two  major  topic  areas:  the  development  of 
the  simulation  model  in  Chapter  III  and  the  illustration  to 
potential  users  of  the  operation  of  the  model  in  Chapter  IV. 
The  development  of  the  system  model  is  traced  within  the 
three  phases  of  the  systems  science  paradigm:  conceptualiza- 
tion,  analysis  and  measurement,  and  computerization.  Both 
structural  and  parametric  models  are  presented  to  aid  in  the 
illustration  of  the  complex  interrelationships  between  sys¬ 
tem  components.  Finally,  the  issues  of  model  verification 
and  validation  are  dealt  with.  During  verification,  the 
model  is  examined  to  ensure  that  it  is  behaving  as  intended. 
This  is  accomplished  by  a  building  block  approach.  First, 
each  part  of  the  model  is  examined  separately  to  verify  its 
correctness.  Then,  the  parts  are  added  one  by  one  until  the 
entire  model  has  been  verified.  During  validation,  a  two¬ 
pronged  approach  is  followed.  The  objective  of  the  first 
phase  of  the  validation  process  is  to  ensure  that  the  model 
coincides  with  the  operation  of  the  real  system.  That  way, 
assumptions  about  the  interrelationships  of  the  model  com¬ 
ponents  are  correct  assumptions  about  the  interrelationships 
of  the  actual  system.  The  objective  of  the  second  phase  is 
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validation  of  the  model  from  an  empirical  point  of  view. 

That  is,  prima  facie  evidence  of  validity  is  gained  if  the 
model  is  used  to  effect  change.  Once  the  model  has  been 
developed,  verified,  and  validated,  the  next  step  is  a  demon¬ 
stration  of  its  operation. 

Chapter  IV  of  this  thesis  contains  an  example  situ¬ 
ation  that  the  model  is  used  to  evaluate.  A  new  RV  is  pro¬ 
posed  and  its  effectiveness  is  evaluated  in  terms  of  its 
success  in  penetrating  a  hypothetical  ballistic  missile 
defense  network  and  destroying  a  target  area.  The  example 
also  contains  the  statistical  experiment  necessary  to  analyze 
such  a  complex  issue.  The  example  is  limited  to  hypothetical 
RVs  and  BMD  systems  because  of  the  security  classification 
of  the  parameters  of  actual  systems.  With  the  User's  Manual 
contained  tr  Appendix  C,  an  analyst  with  access  to  actual 
system  capabilities  could  easily  evaluate  real  systems.  The 
example  simulation  does  demonstrate  some  of  the  capabilities 
of  the  model  and  illustrates  the  complexities  and  significant 
factors  of  the  actual  system.  The  statistical  analysis  leads 
to  the  conclusions  and  recommendations  contained  in  Chapter 
V,  the  final  chapter  of  this  thesis. 


II  The  System 


Introduction 

The  system  under  consideration  in  this  thesis  is  a 
hypothetical  target  area  that  is  defended  by  an  antiballis- 
tic  missile  system  and  attacked  by  ICBM  reentry  vehicles. 

In  order  to  describe  the  system,  the  goals  of  the  system's 
components  are  first  presented.  This  is  followed  by  an 
explanation  of  the  three  major  components  of  the  system,  that 
is,  the  attacking  RVs,  the  defending  ABM  system,  and  the  tar¬ 
get  area.  The  assumptions  made  about  the  components  in  these 
explanations  are  central  to  the  understanding  of  both  the 
actual  system  and  the  model  contained  in  this  thesis. 

The  System  Goals 

Since  all  systems  are  goal  oriented,  the  first  step 
toward  gaining  an  understanding  of  the  system  under  consid¬ 
eration,  and  developing  the  system's  structure,  is  an  iden¬ 
tification  of  the  goals  of  the  system.  There  are  two  oppos¬ 
ing  sides  in  this  system,  the  attacking  RVs,  and  the  defend¬ 
ing  BMD  and  target  area.  Each  side  has  goals.  Figure  3  is 
an  illustration  of  an  assumed  goal  hierarchy  for  the  attack¬ 
ers.  The  intent  of  Figure  3  is  to  express  the  idea  that  the 
attackers  must  first  penetrate  the  defenses,  the  lowest  level 
goal,  and  then  inflict  a  specific  level  of  damage  in  order 
to  accomplish  the  highest  level  goal,  removal  of  the  oppo¬ 
nent's  war-making  capability. 
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defenses  actually  support  the  enemy's  goal  of  offensive 
capability  by  preventing  or  restricting  damage  to  those  offen¬ 
sive  forces.  The  struggle  for  the  attainment  of  these  two 
dichotomous  goal  hierarchies  leads  to  the  structure  of  the 
system  components  as  discussed  in  the  next  sections. 

The  Attackers 

The  attacking  force  is  composed  of  ICBMs  carrying 
actual  reentry  vehicles,  decoys  designed  to  simulate  reentry 
vehicles  or  deceive  BMD  radars,  and  chaff  employed  to  screen 
the  attacking  force  from  enemy  radars.  The  total  system 
composed  of  all  the  ICBMs  attacking  target  areas  similar  to 
the  hypothetical  target  area  can  be  decomposed  into  a  single 
system  made  up  of  the  RVs,  decoys,  and  chaff  carried  on  board 
a  single  ICBM  without  loss  of  meaningful  detail  due  to  the 
explicit  assumption  of  homogeneity.  Each  of  the  four  sub¬ 
systems  of  this  system  will  be  described  in  light  of  their 
functional  contributions  to  the  attainment  of  the  attacking 
force's  goal  hierarchy. 

The  ICBM.  The  United  States  had  a  successful  test 
of  its  first  ICBM,  the  Atlas,  in  the  late  1950s  (Ref  3:75). 
The  first  operational  squadron  was  turned  over  to  the  Strate¬ 
gic  Air  Command  in  September  of  1959.  At  the  same  time,  the 
U.S.  government  conceded  that  the  USSR  had  a  similar  opera¬ 
tional  capability  (Ref  3:112).  The  deployed  systems  of  the 
U.S.  have  evolved  from  those  early  ICBMs  with  single,  large 
yield  warheads  to  the  newest  version  of  the  Minuteman  III 


ICBM  that  has  three,  relatively  low  yield  but  highly  accurate 
RVs  per  ICBM.  The  range  of  strategic  targeting  options  has 
also  evolved,  based  on  the  capabilities  of  the  deployed  sys¬ 
tem.  The  ICBM  of  the  near  future  is  the  MX.  It  is  capable 
of  carrying  many  RVs  with  yields  and  accuracies  sufficient 
to  qualify  them  as  hard  target  killers. 

The  Reentry  Vehicles .  ICBM  reentry  vehicles  can  be 
characterized  as  belonging  to  two  general  categories:  those 
which  have  relatively  poor  accuracy  and  large  yield  making 
them  suitable  for  soft  targets,  and  those  which  have  excel¬ 
lent  accuracy,  low  yield,  and  are  suitable  for  use  against 
hard  targets.  Of  course,  RVs  from  the  second  category  can 
be  used  against  soft  targets  though  this  may  be  inefficient 
from  an  optimal  use  of  weapons  standpoint.  As  a  result,  the 
ICBM  capability  can  be  said  to  be  driving  the  strategy  of 
their  application.  An  example  of  this  Is  the  way  in  which 
the  "mutual  assured  destruction"  strategy  of  the  1960s  gave 
way  to  the  limited  counterforce  proposals  of  Secretary  of 
Defense  Schlesinger  In  the  mid-1970s  (Ref  19:3-4).  This 
shift  in  U.S.  strategic  doctrine  from  targeting  large  cities 
(soft  targets)  to  targeting  enemy  offensive  forces  (hard  and 
soft  targets)  was  possible  because  of  the  development  and 
deployment  of  the  second  category  of  RV  (Ref  19:5-6).  The 
trend  of  RV  research  and  development  has  continued  to  be  in 
the  direction  of  carrying  multiple  RVs  of  the  second,  or 
hard  target  killing,  category.  The  Mk  12 A  reentry  vehicle 


which  is  deployed  on  some  current  ICBMs  is  a  good  example 
of  this  trend,  as  is  the  Advanced  Ballistic  Reentry  Vehicle 
(ABRV)  that  may  be  deployed  in  the  future  (Ref  8;  23).  A 
reentry  vehicle  currently  under  development  may  offer  promise 
as  the  ultimate  in  ICBM-delivered ,  hard  target  killers.  This 
new  RV  is  the  Advanced  Maneuvering  Reentry  Vehicle,  or  more 
simply,  AMARV. 

There  are  two  capabilities  that  make  the  AMARV  dif¬ 
ferent  from  other  RVs.  The  first  is  indicated  by  its  name, 
that  is,  a  maneuvering  capability  after  the  RV  separates 
from  the  post-boost  vehicle.  Unlike  other  RVs  that  fall  bal- 
listically  to  the  target,  the  AMARV,  with  the  use  of  a  digi¬ 
tal  computer  on  board  the  RV,  can  perform  a  series  of  con¬ 
trolled  maneuvers  on  the  way  to  the  target  area.  This  has 
two  purposes.  The  first  is  to  make  it  harder  for  enemy  de¬ 
fenses  to  destroy  the  AMARV.  The  problem  the  BMD  must  solve 
is  predicting  where  an  attacking  RV  is  going  and  then  launch¬ 
ing  a  missile  to  intercept  an  RV  attacking  a  target  that  must 
be  defended.  The  maneuvering  of  the  AMARV  makes  both  the 
trajectory  prediction  and  actual  interception  much  more  diffi 
cult.  This  type  of  AMARV,  termed  an  evader .  compares  the 
position  of  the  AMARV  to  the  prestored  target  coordinates 
and  maneuvers  to  the  coordinates.  The  other  purpose  of  the 
ability  to  maneuver  is  to  support  the  second  capability  of 
the  AMARV  that  makes  it  different  from  the  other  RVs,  termi¬ 
nal  guidance  (Ref  1:43-45;  4). 

The  AMARV  can  theoretically  be  deployed  with  a  target 
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recognition  capability.  This  differs  from  using  prestored 
target  coordinates.  Together  with  the  ability  to  maneuver, 
the  terminally  guided  version  of  the  AMARV  can  theoretically 
be  steered  with  almost  perfect  accuracy  through  the  defenses 
to  the  target.  Both  versions  of  the  AMARV,  and  other  RVs, 
are  designed  to  be  used  in  conjunction  with  penetration  aids, 
such  as  chaff  and  decoys. 

Chaff  and  Decoys.  Chaff  and  decoys  are  designed  to 
aid  the  penetration  of  actual  RVs  through  the  defenses  to  the 
target  area.  Chaff  is  deployed  in  long,  roughly  cylindrical 
tubes.  The  RVs  and  decoys  are  positioned  inside  these  tubes 
and  are  screened  from  the  BMD  radars  by  the  chaff.  The  RVs 
and  decoys  travel  much  faster  than  the  falling  chaff  and, 
therefore,  separate  from  the  chaff  screen.  After  this  sepa¬ 
ration,  another  means  of  deception  is  used. 

One  type  of  decoy,  designed  to  present  the  same 
observable  characteristics  as  actual  RVs  such  as  radar  cross 
section  and  wake,  is  used  to  prevent  the  BMD  radars  from 
correctly  identifying  threats.  This  type  of  decoy  is  called 
a  precision  decoy.  Traffic  decoys  are  smaller  and  are 
deployed  in  large  numbers  in  an  attempt  to  overload  the  BMD 
tracking  and  data  handling  systems  (Ref  7:40).  Thus,  chaff 
and  decoys  are  deployed  to  aid  penetration  by  first  screen¬ 
ing  and  then  disguising  the  actual  RVs. 

The  Defenders 

Ballistic  missile  defense  falls  into  three  basic 


categories,  boost  phase  defense,  midcourse  defense,  and  ter¬ 
minal  defense.  Some  proponents  of  boost  phase  defense  envi¬ 
sion  a  satellite-based,  directed  energy  weapon  that  destroys 
the  ICBM  in  the  period  immediately  after  takeoff  while  it  is 
"large,  slow,  and  soft,"  and  before  multiple  RVs  and  penetra¬ 
tion  aids  can  be  deployed  (Ref  7:^3)*  Boost  phase  defense 
technology  will  probably  not  be  available  in  this  decade  but 
could  become  the  BMD  of  the  1990s  (Ref  7:43-44).  Of  more 
near  term  interest  is  midcourse  defense,  that  is,  defense 
that  takes  place  after  boost  and  before  atmospheric  reentry. 

The  advocates  of  midcourse  defense  call  for  the  use 
of  an  optically  guided  interceptor  that  could  use  nonnuclear 
munitions  as  a  kill  mechanism.  The  BMD  program  is  approach¬ 
ing  the  threshold  of  breaking  the  technology  barrier  of  mid- 
ccurse  defense  according  to  some,  but  deployable,  operational 
systems  are  still  many  years  in  the  future  (Ref  7:41-43). 

Historically,  the  third  type,  terminal  defense,  has 
received  the  most  emphasis  because  the  technology  was  most 
easily  attainable.  Terminal  defense  takes  place  in  the  atmo¬ 
sphere  which  favors  the  defense.  The  atmosphere  slows  down 
attacking  RVs,  filters  out  the  chaff  and  other  penetration 
aids,  and  provides  RV  wake  observables  that  aid  in  threat 
identification  and  discrimination  (Ref  7:38).  All  deployed 
BMD  systems  are  of  this  type  and  all  near  future  systems 
(1980s)  will  probably  be  of  this  type  (Ref  7:42-43).  As  a 
result,  chis  thesis  only  considers  terminal  BMD  systems. 
Terminal  systems  are  composed  of  the  command  and  control 


facilities,  such  as  the  radars  and  computers,  and  the  inter¬ 
ceptors. 


The  Command  and  Control  Facilities .  The  computers 
and  radars  of  a  terminal  defense  network  must  be  extremely 
capable  because  of  the  compressed  time  frame  of  the  engage¬ 
ment.  The  computers  must  be  capable  of  performing  tens  of 
millions  of  instructions  per  second  and  handling  probably 
the  most  demanding  data  processing  requirement  of  any  weapon 
system  (Ref  7:38).  There  are  typically  about  fifteen  seconds 
for  the  system  to  track,  identify,  launch  an  interceptor, 
and  have  time  for  the  interceptor  to  climb  to  the  engagement 
(Ref  7:38).  This  leads  to  the  defense’s  use  of  large,  elec¬ 
tronically  steered,  phased  array  radars  and  extremely  high 
speed  computers  that  are  capable  of  handling  hundreds  of  tar¬ 
gets  at  once.  The  BMD  critic’s  contention  that  these  are 
targets  that  will  be  destroyed  first,  thus  eliminating  the 
defense,  can  be  negated  by  the  concept  of  mobility  that  is 
being  stressed  in  some  current  and  future  BMD  designs.  The 
defense  can  be  mounted  on  vehicles  that  will  prevent  them 
from  being  easily  target able.  The  purpose  of  the  command 
and  control  facilities  is  to  acquire  the  targets,  identify 
the  threats  correctly,  and  launch  interceptors. 

The  Interceptors.  The  interceptors  for  a  terminal 
BMD  are  of  two  varieties,  high  altitude  and  low  altitude. 

The  high  altitude  or  endo-loiter  interceptor  is  launched  when 
the  threat  cloud,  made  up  of  chaff,  RVs,  and  decoys,  is 
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initially  detected,  long  before  individual  threats  are  iden¬ 
tified.  The  endo-interceptor  climbs  to  the  area  where  the 
radar  will  be  able  to  correctly  identify  a  reentry  vehicle 
and  waits  for  that  identification  to  take  place.  It  is  then 
directed  to  the  identified  RV  and  detonated.  Endo-interceptors 
use  relatively  high  yield  weapons  (megaton  range)  and  the 
primary  kill  mechanism  is  X-ray  radiation  throughout  most  of 
the  El  region  because  the  intercepts  take  place  at  high  alti¬ 
tudes  where  air  is  not  dense.  The  X-rays  can  kill  RVs  in  a 
variety  of  ways  such  as  damaging  the  RV  heat  shield  causing 
the  RV  to  either  burn  up  on  reentry  or  become  aerodynamically 
unstable.  The  same  radiations  that  are  used  to  destroy 
attacking  RVs  can  Indirectly  hinder  the  efforts  of  the  defend¬ 
ers  through  a  phenomenon  known  as  blackout. 

Radar  blackout  can  be  caused  by  the  detonation  of  the 
ABM  warhead.  That  much  is  certain.  The  problem  is  compounded, 
however,  by  other  uncertainties.  There  is  no  clear  analytic 
approach  to  the  problem  and  the  bans  on  atmospheric  testing 
preclude  experimentation.  As  a  result,  no  one  knows  exactly 
what  the  effects  are  or  how  much  degradation  the  BMD  system 
effectiveness  will  suffer.  Any  effects  will  be  magnified 
because  of  the  large  yield  weapons  used  in  the  endo-intercept 
region  (Ref  M). 

The  low  altitude  or  point  interceptor  uses  a  smaller 
weapon  (several  kiloton  range)  and  is  launched  to  destroy 
any  RVs  that  evade  the  endo-interceptors  or  after  the  inven¬ 
tory  of  endo-interceptors  has  been  exhausted.  Because  X-rays 
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are  more  readily  attenuated  by  the  denser  atmosphere  of  the 
terminal  region  than  neutrons,  the  primary  kill  mechanism  in 
the  point  intercept  region  is  neutron  radiation.  Neutrons 
destroy  RVs  by  causing  damage  to  the  Internal  electronic 
equipment  or  to  the  RV  warhead  itself.  Those  reentry  vehicles 
that  evade  the  defense  continue  to  the  target  area. 

The  Target  Area 

The  hypothetical  target  area  for  this  thesis  is 
assumed  to  be  made  up  of  targets  requiring  several  thousands 
of  pounds  per  square  inch  (psi)  of  overpressure  to  be  de¬ 
stroyed,  or  so-called  hard  targets.  A  good  example  of  this 
kind  of  target  area  is  an  ICBM  squadron  with  its  command  and 
control  facilities  and  ICBMs  in  hardened  silos. 

Summary 

The  system  under  observation  is  composed  of  the 
attacking  RVs,  with  the  associated  chaff  and  decoys ;  the 
defending  BMD  system,  made  up  of  endo-im;erceptors  and  point 
interceptors;  and  the  target  area.  Because  the  actions  of 
any  system  are  very  often  different  than  the  sura  of  the  ac¬ 
tions  of  the  parts,  nothing  is  gained  by  studying  only  the 
parts.  The  interrelationships  between  the  parts  must  also 
be  studied.  This  study  can  be  accomplished  with  a  simulation 
model.  The  model  discussed  in  the  next  chapter  will  be  used 
to  analyze  both  the  interactions  of  the  system  components  and 
the  components  themselves. 
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Ill  The  Model 


Introduction 

This  chapter  presents  the  translation  process  from 
the  actual  system  discussed  earlier  to  the  simulation  model 
of  the  system.  The  development  of  the  model  follows  the  sys¬ 
tem  science  paradigm  of  conceptualization,  analysis  and  mea¬ 
surement,  and  computerization.  In  the  conceptualization 
phase,  the  components  of  the  system  are  identified  and  the 
first  attempts  at  identifying  their  interrelationships  are 
made.  The  central  goal  of  this  phase  is  a  structural  model 
of  the  system.  To  accomplish  this,  the  attackers,  defenders, 
and  the  target  area  are  closely  studied  with  the  objective 
of  gaining  a  complete  understanding  of  how  the  actual  system 
functions.  This  understanding  of  the  system  is  vital  for 
the  second  phase  of  model  development,  analysis  and  measure¬ 
ment  . 

The  analysis  and  measurement  phase  begins  the  quanti¬ 
fication  of  the  model.  The  underlying  mathematics  and  physics 
of  the  engagements  between  the  attacking  RVs  and  the  endo- 
interceptors  and  point  interceptors  and  the  damage  to  the 
target  area  by  penetrating  RVs  are  presented.  These  form 
the  parametric  model  of  the  system  that  is  necessary  for  the 
third  phase  of  the  paradigm,  computerization. 

The  objective  of  the  computerization  phase  of  model 
development  is  the  combination  of  the  structural  and  para¬ 
metric  models  to  form  a  computer  simulation  model  of  the  actual 
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system.  The  discussion  of  the  computerization  phase  includes 
the  main  program  and  supporting  subroutines.  The  detailed 
and  annotated  codes  are  contained  in  Appendix  C. 


The  final  phase  of  the  model  development  is  the  veri¬ 
fication  and  validation  of  the  model.  The  last  section  of 
this  chapter  details  the  efforts  designed  to  ensure  that  the 
model  behaves  as  intended  and  reflects  a  realistic  view  of 
the  real  system. 

Conceptualization 

The  conceptualization  phase  of  model  development 
begins  with  the  definition  of  the  goals  of  the  simulation. 
These  are  outlined  in  the  first  chapter  of  this  thesis.  The 
major  objective  is  to  provide  a  model  that  permits  analysis 
in  all  areas  of  the  engagements  between  the  attacking  RVs 
and  the  defending  ABM  system.  In  order  to  accomplish  this, 
the  actual  system  must  be  studied  and  the  components  and 
functional  relationships  identified.  The  three  major  com¬ 
ponents  and  their  capabilities  are  defined  as  precisely  as 
possible  so  that  they  can  be  modeled.  The  next  sections 
detail  the  thesis  definitions  of  and  assumptions  concerning 
the  major  components. 

The  Attackers.  The  three  subsystems  that  comprise 
the  attacking  force  are  the  ICBMs,  the  RVs,  and  the  penetra¬ 
tion  aids  (chaff  and  decoys).  The  overall  battle  can  be 
modeled  on  the  level  of  one  ICBM  and  its  load  of  RVs  and 
penetration  aids  versus  one  BMD  system.  The  hypothetical 


ICBM  of  this  thesis  can  carry  up  to  4000  pounds  of  RVs  and 
penetration  aids.  It  is  a  solid  fuel  ICBM  and  has  a  range 
with  the  defined  mission  weight  comparable  to  the  newest 
U.S.  systems(Ref  26).  This  thesis  does  not  model  the  launch 
of  the  ICBM.  Such  issues  as  prelaunch  survivability  and 
launch  reliability  are  not  addressed.  Rather,  the  simulation 
begins  as  the  post-boost  vehicle  has  reached  the  point  at 
reentry  when  the  RVs  and  aids  will  be  deployed.  The  post¬ 
boost  vehicle  deploys  the  RVs,  one  by  one,  with  some  time 
between  separation  of  each  RV.  This  time  between  RV  deploy¬ 
ments  is  affected  by  factors  such  as  the  desired  spacing 
between  objects  in  a  chaff  cloud,  and  the  ground  distance 
between  targets. 

The  RVs  that  can  be  carried  by  the  hypothetical 
booster  can  range  from  RVs  such  as  are  presently  deployed 
on  Minuteman  III,  to  those  planned  for  deployment  on  a  future 
MX  ICBM,  to  theoretical,  as  yet  untested  RVs.  For  demonstra¬ 
tion  purposes  in  the  model,  this  thesis  uses  two  different 
types,  a  maneuvering  RV  (MRV)  and  a  nonmaneuvering,  or  bal¬ 
listic  RV  (BRV).  The  assumed  capabilities  are  displayed  in 
Table  I. 


TABLE  I 

Model  Reentry  Vehicle  Performar  ce  Characteristics 

RV  Type  Yleld(kt)  CEP (feet)  Weight (lbs) 

MRV  150  200-600 

600 


BRV 


350 
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350 

350 


The  values  in  Table  I  are  assumed  for  demonstration 
purposes  only  and  do  not  reflect  the  performance  character¬ 
istics  for  any  actual  weapon  system.  The  values  can  be 
changed  with  little  difficulty  to  allow  the  analysis  of  an 
actual  system.  A  demonstration  of  the  ease  of  changing 
parameters  is  found  in  the  example  experiment  found  in  Chap¬ 
ter  IV. 

Both  of  these  hypothetical  RVs  are  assumed  to  be  pro¬ 
tected  against  nuclear  effects  from  interceptor  weapon  deto¬ 
nations  to  some  designed  vulnerability  levels.  These  levels 
are  changeable  in  the  model.  The  RVs  are  electrically  com¬ 
patible  and  similar  in  size  so  that  they  can  be  placed  at  any 
RV  station  in  the  post-boost  vehicle.  This  means  that  the 
two  types  are  interchangeable  but  the  maximum  total  number 
of  RVs  that  can  be  carried  at  one  time  is  ten,  regardless  oi 
whether  or  not  penetration  aids  are  deployed. 

Although  the  actual  deployment  of  penetration  aids 
is  not  simulated  in  the  model,  their  use  is  an  implicit 
assumption.  This  assumption  results  in  the  150,000  feet 
altitude  that  represents  the  level  at  which  individual  RVs 
are  correctly  identified  by  the  ABM  radars.  Without  the  use 
of  penetration  aids,  the  discrimination  altitude  would  be 
much  higher,  resulting  in  a  much  larger  engagement  regime. 

The  consideration  of  either  traffic  decoys  designed  to  over¬ 
load  BMD  command  and  control  subsystems  or  decoys  designed 
to  give  the  same  observable  characteristics  as  actual  RVs  is 
also  reflected  in  the  probabilities  of  correct  discrimination 
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in  the  model  networks  for  the  El  and  point  defenses.  The  use 
of  electronic  countermeasure  (ECM)  decoys  designed  to  jam  or 
deceive  BMD  radars  can  also  be  modeled  by  reflecting  the  im¬ 
pact  of  their  use  on  the  probability  of  correctly  discrimi¬ 
nating  the  threats  by  the  BMD. 

Just  as  the  attacking  forces  can  be  represented  by 
the  subsystems  of  a  single  ICBM  and  its  load  of  RVs  and 
penetration  aids,  the  defender^  can  be  represented  by  a  sin¬ 
gle  system  made  up  of  command  and  control,  endo-interceptor 
and  point  interceptor  subsystems. 

The  Defenders.  The  performance  of  the  command  and 
control  subsystem  of  the  BMD  is  modeled  in  terms  of  its  cap¬ 
ability  to  correctly  discriminate  and  engage  the  threat. 

This  performance  forms  a  part  of  the  probabilities  of  dis¬ 
crimination  used  in  the  model  network  to  route  entities  to 
engagements  if  discriminated  correctly  or  to  bypass  the 
engagement  simulation  if  not  correctly  discriminated.  How¬ 
ever,  the  probabilities  of  discrimination  figures  reflect 
an  aggregation  of  many  more  probabilities.  For  example,  the 
launch  reliabilities  of  the  interceptors  are  included  as  well 
as  the  probabilities  of  weapon  fusing  and  detonation. ^  As  a 
result,  the  probabilities  of  discrimination  assessed  for  the 
BMD  system  are  critical  and  the  sensitivity  of  the  system  to 
these  values  is  analyzed.  However,  the  actual  calculation' 
of  these  figures  is  outside  the  scope  of  this  thesis.  This 
same  type  of  simplifying  assumption  is  used  in  the  simulations 


of  the  er.gag  nents  between  the  interceptors  and  the  attack¬ 
ing  reentry  vehicles. 

The  underlying  aerodynamic  performance  of  the  inter¬ 
ceptors  is  not  modeled.  The  concept  of  weapon  radius  (WR) 
is  used  in  the  simulation  of  engagements.  The  mathematical 
derivation  of  the  weapon  radius  is  contained  in  Appendix  B 
but,  in  general  terms,  the  weapon  radius  is  used  to  generate 
the  probabilities  3f  kill  in  the  engagements.  The  weapon 
radius,  of  course,  varies  based  on  the  specific  interceptor, 
endo-interceptor  or  point  interceptor. 

The  hypothetical  endo-interceptor  is  assumed  to  have 
some  form  of  loiter  capa~~lity.  It  is  launched  well  in  ad¬ 
vance  of  actual  target  discrimination  and  climbs  to  the  upper 
level  of  the  endo-incerr  -pt  region.  Once  the  final  discrimi¬ 
nation  has  taken  place,  ;he  .^ao-interceptor  is  directed  to 
the  engagement  with  the  attacking  thrQat.  Since  this  engage¬ 
ment  takes  place  in  atmospheric  regions  of  ^ow  air  density, 
the  powered  maneuvering  of  the  endo-interceptor  can  partially 
negate  the  aerodynamic  maneuvering  of  the  MRV.  Rather  than 
model  the  intricate  maneuvering  process  engendered  by  the 
engagement  process,  the  model  uses  the  weapon  radius  concept. 
The  weapon  radius  is  affected  by  a  variety  of  factors  that 
are  illustrated  in  Appendix  B.  The  warhead  on  the  inter¬ 
ceptor  is  a  large  yield,  thermonuclear  device  in  the  megaton 
range.  The  intent  of  the  use  of  this  large  warhead  is  to 
overcome  the  disadvantages  of  the  interceptor  in  the  endo- 
intercept  engagement.  One  of  these  disadvantages  is  the 
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classical  one  of  the  defense  being  forced  to  react  to  the 
initiative  of  the  attacker.  Since  the  attacker  is  still  far 
above  the  target,  an  attacking  maneuvering  RV  can  perform 
violent  maneuvers  in  the  endo-intercept  region  and  still  have 
enough  time  to  maneuver  to  the  target.  This  is  less  true  in 
the  point  intercept  region. 

The  interceptor  in  the  point  intercept  region  is 
assum  ;d  to  be  a  small,  high  acceleration,  extremely  maneu¬ 
verable,  solid  fuel  missile  with  a  warhead  in  the  several 
kiloton  range.  This  smaller  weapon  is  needed  partly  because 
of  less  time  to  maneuver  for  an  attacking  RV  and  partly  to 
prevent  collateral  damage  from  the  detonation  of  the  ABM.,  A 
related  issue  to  this  collateral  damage  problem  is  the  poten¬ 
tial  for  ABM  radar  blackout  as  a  result  of  weapon  detonation. 
This  is  not  modeled  in  this  thesis  because  no  concrete  infor¬ 
mation  dealing  with  the  blackout  potential  could  be  found. 

In  addition  to  the  assumed  capabilities  of  the  hypo¬ 
thetical  interceptors,  the  strategy  of  the  defense  is  assumed 
to  be  to  engage  each  RV  with  one  ABM  at  a  time.  This  means 
there  is  no  netting,  or  engaging  a  single  RV  with  two  or  more 
interceptors  from  a  single  zone  (endo  or  point).  However, 
a  reentry  vehicle  that  successfully  penetrates  the  endo- 
intercept  zone  can  be  engaged  by  a  point  interceptor.  An 
additional  assumption  is  that  an  interceptor  can  only  damage 
the  specific  RV  it  is  engaging;  collateral  damage  of  a  RV 
other  than  the  one  targeted  is  not  modeled. 


33 


The  Target  Area.  The  target  area  used  in  this  model 
is  made  up  of  fixed,  homogeneous,  hard  targets.  The  model 
contains  the  explicit  assumption  that  these  targets  are  spaced 
so  that  each  individual  target  can  only  be  damaged  by  the 
specific  RV  attacking  it.  That  is,  one  RV  missing  its  target 
cannot  damage  another  target.  An  example  of  this  type  of 
target  area  is  an  ICBM  squadron. 

The  FORTRAN  subroutines  that  calculate  the  level  of 
damage  suffered  by  each  target  only  consider  overpressure. 
There  are  two  reasons  for  this.  The  first  is  ease  of  calcu¬ 
lation.  Although  ICBM  silos,  for  example,  are  susceptible 
to  damage  caused  by  ground  shock,  the  calculation  of  the  lev¬ 
el  of  damage  caused  by  ground  shock  is  an  imprecise  science 
at  best.  The  second  reason  for  only  considering  overpressure 
effects  is  that  damage  caused  by  these  effects  can  be  visu¬ 
ally  confirmed  by  reconnaissance  after  the  attack. 

Another  assumption  in  the  model  is  that  each  target 
i3  attacked  by  only  one  RV.  There  is  no  N-on-one  strategy 
which  would  have  two  or  more  RVs  attacxing  a  single  target. 

The  three  major  components  discussed  previously,  the  attackers 
the  defenders,  and  the  target  area,  form  the  basis  for  the 
structural  model  of  the  system. 

The  Structural  Model.  The  development  of  a  struc¬ 
tural  model  of  the  system  being  studied  is  vital  in  order 
to  gain  a  fuller  understanding.  The  three  components  already 
identified,  the  attackers,  the  defenders,  and  the  target  area. 


can  be  used  to  develop  a  basic  causal  loop  diagram,  the  begin¬ 
ning  of  structural  model  development. 

The  intent  of  Figure  5  is  to  visually  demonstrate 
the  interactions  between  components.  For  example,  the  arrow 
between  one  defenders  and  the  target  area  indicates  that  a 
large  number  of  defenders  should  result  in  a  negative  effect 
on  the  target  area  component  as  measured  by  the  target  damage 
level,  that  is,  less  damage.  Similarly,  an  increase  in  the 
number  of  attackers  should  result  in  a  positive  effect  or 
more  target  area  damage.  This  simplified  view  of  the  system 
can  be  made  more  complex  with  an  introduction  of  some  of  the 
parameters  involved.  Table  II  is  a  partial  listing  of  the 
system  parameters. 

Although  these  parameters  do  complicate  the  issue, 
the  actual  structure  of  the  system  is  still  relatively  straight¬ 
forward.  The  underlying  structure,  or  logic  flow,  of  the  sys¬ 
tem  is  illustrated  in  Figure  6.  The  figure  uses  the  four 
parts  of  the  system  as  described  by  Matlin,  the  weapon  com¬ 
plex,  the  target  complex,  the  engagement  model,  and  the  dam¬ 
age  model  (Ref  12:337). 
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TABLE  II 


The  Missile  Allocation  Problem:  A  Partial 
Parameter  Listing  (Ref  12:335) 


ATTACKER 

(for  each  weapon  type) 
Availability 
Survivability 


Booster  Characteristics 
Location 
Range 
Accuracy 
Memory 
Penetration  aid  mix  and 
effectiveness 

RV  Characteristics 
Accuracy 
Yield 

A  range,  A  crossrange 
capability 

Command  and  control 
Target  sharing 
Reprogramming 
Weapon  commitment 


Reliability 

Bomb  damage  assessment 

Etc . 


Signature 

Indirect  bomb  damage  assessment 

Survivability 

Etc . 


Redundancy  and  survivability 

Reaction  time 

Etc. 


OBJECTIVES 


Offense 

Damage  criterion  per  target 
Total  criteria 

Defense 

Damage  criterion  per  target 
Over-all  criteria 


INTELLIGENCE 

What  the  defense  knows  about  the  offense 
Number  of  attacking  missiles 
Number  of  salvos 
Salvo  composition 
Missile  payload  mix 
Etc. 

What  the  offense  knows  about  the  defense 

Terminal  ABM  inventories,  distribution  among  targets 
Area  ABM  inventories 
Preferential  defense  strategy 
Etc, 
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TABLE  II  -  Continued 


DEFENDER 


Target  complex 

Locations/uncertainties 

Configurations 

Values 

Hardness 

Defenses 

Terminal  inventories 
Terminal  effectiveness 
Area  coverages 

Area  characteristics  (inventory,  choice,  effectiveness) 

Radar  capabilities  (tracking,  impact  prediction,  capability) 
Civilian  defense  program 

SCENARIO  (strategic  environment) 

Enemy  first  strike 

Assured  destruction  mission  for  U.S. 

U.S.  attrition 

U.S.  first  strike 

Damage  limiting  mission  for  U.S. 

(Counterforce,  countervalue,  retained) 
mix 

Selective  threat  targeting  philosophy 
Progressive  confrontation  (tit-for-tat )  targeting 
Collateral  damage  restrictions 
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Pig  6.  Underlying  Logic  Plows  (Ref  12:336) 


Using  the  aforementioned  list  of  system  parameters 
and  the  underlying  logic  flow  of  Matlin,  the  basic  causal 
loop  diagram  can  be  expanded.  Combining  the  knowledge  of 
the  system  components  and  their  interactions  gained  from  the 
basic  causal  loop  diagram,  the  list  of  system  parameters, 
and  Matlin' s  underlying  logic  flow  yields  the  structural 
model  of  the  system.  This  is  depicted  in  Figure  7. 

An  integral  part  of  the  structural  model  is  an  iden¬ 
tification  of  the  system  variables  by  type.  The  control 
variables,  stochastic  variables,  and  response  variables  must 
be  Identified.  This  identification  is  depicted  in  Table  III. 
The  control  variables  are  those  input  variables  over  which 
the  decision  maker  has  control.  This  control  may  exist 
through  policy  decisions  such  as  equipment  design  and  devel¬ 
opment  decisions,  and/or  employment  decisions.  Depending 
upon  the  purpose  for  which  the  model  is  used,  the  decision 
maker  may  have  to  rely  on  intelligence  and  technological 
estimates  to  control  these  variables.  Stochastic  variables 
are  those  input  variables  over  which  the  decision  maker  has 
no  control.  The  values  of  these  variables  are  probabilistic. 
The  response  variable  is  the  output  variable  that  is  used  to 
measure-  the  effectiveness  of  a  component  or  the  system. 

The  structural  model  illustrates  the  framework  of  the 
system.  The  components  of  the  system  are  apparent,  as  are 
their  interrelationships.  In  order  to  measure  these  inter¬ 
actions,  however,  the  model  must  be  quantified.  Some  numeri¬ 
cal  relationships  must  be  developed.  These  mathematical 
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relationships  form  a  parametric  model,  the  central  goal  of 
the  next  phase  of  model  development,  analysis  and  measure¬ 
ment  . 

Analysis  and  Measurement 

The  understanding  of  the  system  developed  in  the  con¬ 
ceptualization  phase  of  system  analysis  is  vital  for  the  de¬ 
velopment  of  the  parametric  model  that  is  the  objective  of 
this  phase. 

The  Parametric  Model.  The  mathematical  relationships 
of  interest  fall  in  four  areas:  (1)  the  time/velocity  func¬ 
tions  for  RV  reentry,  (2)  the  X-ray  and  neutron  calculations, 
(3)  the  weapon  radius  (WR)  calculation*1  used  to  determine 
probabilities  of  RV  kill  by  the  defense,  and  (4)  the  target 
damage  level  calculations.  The  last  three  are  contained  in 
Appendices  A  and  B.  There  are  several  assumptions  central 
to  the  development  of  the  parametric  model.  These  are  made 
explicit  as  they  are  needed.  Figure  8  is  used  to  develop 
the  time/velocity  functions  for  RV  reentry.  If  a  straight 
line  RV  trajectory  is  assumed,  that  is  no  maneuvering,  then 
X  represents  the  RV  reentry  angle.  Since, 


then. 


sin  (X) 


vertical  distance  through  SI  zone 
slant  distance  traversed  by  RV 
through  El  zone 


(1) 


the  slant  distance  traveled  80,000 

by  the  RV  through  the  El  zone  sin(X) 
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Pig  8.  Sample  RV  Trajectory 


This  assumes,  of  course,  that  there  are  exact  upper  and  lower 
boundaries  for  the  El  zone.  The  same  assumption  applies  for 
the  PI  zone  and,  as  a  result, 

the  slant  distance  traveled  M  60,000 
by  the  RV  through  the  PI  zone  sin(X) 

Eqs  (2)  and  (3)  lead  to  the  calculation  of  time  spent  in  each 
region  with  the  assumption  of  constant  velocity,  V. 

Time  to  traverse  El  zone  *  .^ant,  Distance  (4) 

Time  to  traverse  PI  zone  =  ----- (5) 
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The  absolute  minimum  time  that  can  elapse  between  the 
time  a  reentry  vehicle  enters  a  zone  and  is  engaged  is  the 
minimum  reaction  time  of  the  interceptor.  The  maximum  time 
from  entry  to  engagement  is  the  total  time  spent  in  the  zone. 
Since  the  hypothetical  ABM  system  of  this  model  is  not  an 
actual  system,  no  data  exist  on  the  probability  distribution 
function  that  describes  the  distribution  of  intercepts  within 
a  zone.  Although  the  natural  tactic  of  the  defense  would  be 
to  engage  as  high  in  each  zone  as  possible,  thereby  maximiz¬ 
ing  the  lethal  radius  because  of  lower  air  density,  the  de¬ 
fense  must  react  to  the  offense.  This  means,  for  example, 
that  in  the  El  zone,  the  ABM  must  accelerate  and  then  maneuver 
to  engage  the  RV.  As  a  result,  relatively  few  engagements 
should  occur  at  the  upper  boundary.  Since  the  RVs  are  capable 
of  much  higher  velocities  than  the  ABMs,  relatively  few  engage¬ 
ments  should  occur  at  the  lower  boundary  as  well  because  this 
implies  a  tail-chase  not  winnable  by  the  interceptor.  The 
set  of  candidate  probability  distribution  functions  to  des¬ 
cribe  the  distribution  of  intercept  altitudes  throughout  each 
region  is  large.  It  includes  a  triangular  form;  a  lognormal; 
a  normal,  either  symmetric  or  skewed;  or  a  beta  probability 
distribution  function,  to  name  just  a  few.  Historical  data 
are  scarce  because  of  the  short  history  of  deployed  BMD  sys¬ 
tems.  Therefore,  analysis  of  the  results  of  specific  inter¬ 
cept  distributions  was  required.  A  test  was  designed  using 
the  weapon  system  capabilities  and  vulnerabilities  for  both 
the  offense  and  defense  as  contained  in  the  example  simulation 


in  Appendix  C.  That  is,  the  vulnerabilities  of  the  attacking 
RV  to  X-ray  and  neutron  fluence  were  set  at  the  same  levels 
as  the  hypothetical  Maneuvering  RV  (MRV).  These  were  70  cal/ 
cm2  and  100  cal/cm2  for  X-ray  sure  safe  and  sure  kill  levels, 
respectively,  and  lx  10 neutrons/cm2  and  1  x  lO1^  neutrons/ 
cm2,  for  neutron  fluence.  The  two  hypothetical  interceptors 
had  yields  of  one  megaton  for  the  endo-interceptor  and  twenty 
kilotons  for  the  point  interceptor.  The  CEP  for  the  endo- 
interceptor  was  3000  feet;  the  CEP  for  the  point  interceptor 
was  1000  feet. 

The  test  consisted  of  forty  engagements  between  a  RV 
and  an  interceptor,  twenty  in  the  El  region  and  twenty  in 
the  PI  region.  The  RV  was  assumed  to  be  correctly  discrimi¬ 
nated.  A  triangular  function  was  used  to  distribute  the 
intercept  altitudes  with  the  most  likely  altitude  set  at  a 
point  most  disadvantageous  to  the  defense.  Because  of  atmo¬ 
spheric  attenuation,  the  lethal  radii  of  both  primary  kill 
mechanisms.  X-ray  and  neutron  fluence,  shrink  as  altitude  de¬ 
creases  and  the  density  of  the  atmosphere  increases.  There¬ 
fore,  the  most  disadvantageous  altitudes  for  the  defense  are 
the  lowest  in  each  zone.  Even  with  the  defense  at  this  disad¬ 
vantage,  all  forty  engagements  resulted  in  the  RV  being  killed. 

As  a  result  of  this  test,  and  other  preliminary  anal¬ 
ysis,  a  normal  probability  distribution  function  was  assumed 
for  the  distribution  of  intercept  altitudes  in  each  zone. 

The  normal  form  is  assumed  to  have  the  following  parameters. 


The  mean, 

(time  to  fly  zone  -  ABM  minimum  reaction  time  (6) 
and  standard  deviation. 


The  probability  distribution  function  can  be  changed  within 
the  model  to  reflect  differences  in  BMD  systems.  Other  func¬ 
tional  forms  can  be  substituted  if  desired. 

Eq  (8)  follows  from  Eqs  (1)  -  (5). 

Time  to  fly  Dead  Zone  *  ~gi nX't)'  (8) 

As  a  result  of  the  calculation  of  the  time  of  inter¬ 
cept  in  each  zone,  the  Intercept  altitudes  can  be  determined. 
The  results  of  Eqs  (9)  and  (10)  are  needed  to  determine  radia^ 
tlon  levels  and  probabilities  of  kill. 

El  altitude  =  150,000  -  (intercept  time) (v) (sin(x) )  (9) 

PI  altitude  »  70,000  -  (intercept  time) (v) (sin(x) )  (10) 

Computerization 

The  objective  of  the  third  phase  of  system  analysis, 
computerization,  is  a  computer  model  of  the  system.  The 
structural  model  developed  during  the  conceptualization  phase 
can  be  changed  into  a  structural  model  more  suitable  for  com¬ 
puterization  by  adopting  a  network  orientation.  By  this,  it 
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is  meant  that  each  RV  must  follow  a  certain  path  from  the 
drop-off  point  from  the  post-boost  vehicle,  through  the  endo- 
intercept  and  point  intercept  regions  and  the  dead  zone,  to 
the  target  area.  Along  the  path  are  decision  points,  or 
nodes,  where  decisions  must  be  made  as  to  whether  or  not  the 
RV  will  continue.  For  example,  in  the  endo-intercept  region, 
the  RV  may  be  engaged  by  the  interceptor.  This  engagement 
can  be  represented  by  a  node  in  a  network  where  the  decision 
is  made  on  the  success  of  the  engagement.  Figure  9  repre¬ 
sents  the  computer  logic  of  this  network. 

The  Simulation  and  Modeling  Language  (SLAM)  was  cho¬ 
sen  for  the  main  program  because  of  its  amenability  to  pro¬ 
gramming  network  structures.  It  provides  a  structural  model 
of  the  system,  laid  out  as  a  network,  and  permits  the  very 
real  advantage  of  control  over  this  structural  model  during 
the  simulation.  This  means  that  not  only  can  parameters  such 
as  weapon  yields  and  accuracies  be  altered,  but  the  physical 
structure  of  the  system  can  be  changed,  all  without  inter¬ 
rupting  the  operation  of  the  model.  The  reason  for  this  tre¬ 
mendous  flexibility  is  the  SLAM  global  variable.  Although 
more  detail  is  presented  in  Appendix  C,  briefly,  global  vari¬ 
ables  car.  be  used  to  open  and  close  network  paths,  activate 
and  deactivate  defenses,  or  redefine  weapon  system  capabil¬ 
ities  and  vulnerabilities,  while  continuing  the  simulation. 
The  use  of  this  special  purpose  language,  and  the  resulting 
flexibility  in  the  simulation  of  a  critically  important 
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Fig  9.  Computer  Model  Flowchart 


subject  area,  is  the  key  to  the  value  of  this  RV/ABM  engage¬ 
ment  model. 

The  language  of  the  subroutines  that  support  the  SLAM 
main  program  is  FORTRAN  V.  These  subroutines  are  the  comput¬ 
erization  of  the  parametric  models  developed  in  the  analysis 
and  measurement  phase.  They  provide  RV  reentry  time  and 
velocity  information,  as  well  as  the  results  of  the  RV/ABM 
engagement  simulations  and  target  damage  levels.  Annotated 
listings  and  detailed  explanations  of  the  main  and  supporting 
programs  are  contained  in  Appendix  C. 

The  end  result  of  the  computerization  phase  is  a  com¬ 
puter  model  that  combines  the  structural  model  of  the  con¬ 
ceptualization  phase  and  the  parametric  model  developed  in 
the  analysis  and  measurement  phase.  Hie  model  contained  in 
this  thesis  allows  detailed  analysis  of  the  system  and  iden¬ 
tification  of  statistically  critical  components,  interactions, 
and  parameters.  But  the  mere  development  of  a  simulation 
model  is  not  enough.  The  model  must  be  verified  and  validated 
before  use  as  an  analytical  tool. 

Verification  and  Validation 

Verification  is  the  process  of  ensuring  that  the  parts 
of  the  model  and  the  entire  model,  itself,  operate  as  intended 
b„  the  designers.  Validation,  on  the  other  hand,  is  ensur¬ 
ing  that  the  model  is  an  accurate  representation  of  the  actual 
system.  Both  processes  must  be  accomplished  if  the  model  is 
to  be  used  to  analyze  the  real  system. 
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The  Verification.  The  verification  process  for  this 
model  was  conducted  in  two  stages,  one  for  the  main  program 
and  one  for  the  supporting  user  functions  (subroutines).  The 
verification  of  ti.e  main  network  was  accomplished  with  a 
building  block  approach.  Each  part  of  the  model  was  run 
alone  with  the  SLAM  TRACE  option  to  make  sure  that  the  model 
part  was  behaving  as  intended.  By  this,  it  is  meant  that 
each  entity  flowed  through  the  network,  gained  attribute  val¬ 
ues  where  and  when  intended,  and  followed  the  designed  path. 
Figure  10  is  an  example  of  the  verification  of  the  first 
block  in  the  network.  With  the  TRACE  option  activated,  all 
entities  (RVs)  could  be  tracked  through  the  Arriving  RVs  block 
of  the  model.  Attributes  3  through  6,  assigned  by  subrou¬ 
tines,  could  be  verified,  and  the  designers  of  the  model 
could  ensure  that  the  block  was  operating  as  intended.  After 
the  first  block  had  been  verified,  the  second,  or  Endo  De¬ 
fense,  block  could  be  verified  separately  and  then  added  to 
the  first  block  and  verified  jointly.  This  process  contin¬ 
ued  until  each  block  had  been  separately  verified  and  the 
main  program  had  been  verified  as  a  whole. 

The  second  stage  of  the  verification  process  involved 
the  supporting  FORTRAN  code.  These  subroutines  supply  infor¬ 
mation  to  be  used  in  the  main  program  such  as  neutron  and 
X-ray  fluence,  or  radiation  levels,  and  damage  levels.  Ana¬ 
lytic,  closed  form  solutions  exist  for  these  calculations 
and  form  the  basis  for  the  programming  in  the  subroutines. 

The  subroutines  were  run  with  test  data  and  their  outputs 
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Arriving  RVs  Verification 


compared  with  the  results  of  the  analytical  solutions.  An 
example  of  this  is  in  Appendix  D.  The  results  of  both  meth¬ 
ods  agreed.  The  accuracy  of  the  analytical  solutions  them¬ 
selves,  and  the  resulting  subroutines,  was  further  verified 
as  a  result  of  their  examination  by  a  recognized  expert  in 
the  field  of  nuclear  weapons  effects,  Dr.  Charles  J.  Bridgman 
of  the  Department  of  Nuclear  Engineering,  Air  Force  Institute 
of  Technology.  Dr.  Bridgman  is  an  expert  in  both  the  theo¬ 
retical  side  of  the  weapons  effects  field  and  the  operational 
side,  weapons  testing,  as  demonstrated  by  his  academic  cre¬ 
dentials  and  his  years  of  experience  in  the  field,  actually 
testing  nuclear  weapons.  Dr.  Bridgman’s  examination  was 
invaluable  not  only  in  mod^l  verification  but  also  in  valida¬ 
tion. 

The  Validation.  The  process  of  model  validation,  or 
ensuring  the  model  accurately  represents  the  actual  system, 
was  based  on  two  definitions  of  what  constitutes  a  valid 
model.  The  first,  and  more  conventional,  definition  states 
that  a  model  is  valid  if  it  accurately  predicts  the  condi¬ 
tions  of  the  actual  system.  Based  on  the  verification  pro¬ 
cess  previously  described,  during  which  an  expert  confirmed 
that  the  output  of  the  model  agreed  with  the  results  reached 
by  widely  accepted,  analytic  solutions,  the  authors  of  this 
thesis  believe  that  the  model  is  valid  as  they  have  defined 
model  validity  in  the  first  sense. 

The  second  definition  of  model  validity  corresponds 
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to  the  position  espoused  by  Naylor,  et.  al.,  in  Computer  Sim¬ 
ulation  Techniques ,  wh&.t  validating  policy  models.  They  feel 
that  since  the  objective  of  a  policy  model  is  to  recommend  a 
specific  course  of  action,  or  policy,  the  adoption  of  that 
policy  by  the  decision-making  authority  validates  the  model 
(Ref  16:315).  Put  more  simply,  if  the  objective  of  a  model 
is  analysis  that  results  in  change,  and  change  results,  the 
model  is  valid.  A  related  position  may  be  offered.  If  the 
objective  of  a  model  is  to  provide  a  tool  for  analysis,  and 
experts  in  the  field  are  willing  to  use  the  tool  (or  model), 
then  the  model  is  valid.  According  to  this  empirical  defini¬ 
tion,  the  model  in  this  thesis  is  valid.  The  assumptions, 
processes,  programming,  and  methods  used  in  this  model  were 
presented  to  USAP  officers  involved  in  both  the  fields  of  ABM 
intelligence  analysis  and  ABM/RV  system  analysis.  They  con¬ 
curred  with  the  model  and  expressed  interest  in  adopting  it 
as  an  analysis  technique  for  the  Strategic  Air  Command  (Ref  22). 

Summary 

This  chapter  has  traced  the  development  of  the  model 
through  the  three  phases  of  the  system  science  paradigm,  that 
is,  conceptualization,  analysis  and  measurement,  and  comput¬ 
erization.  It  also  detailed  the  efforts  of  the  authors  to 
both  verify  and  validate  the  model.  The  next  chapter.  Chapter 
IV,  contains  an  example  problem  designed  to  demonstrate  the 
operation  of  the  model  and  permits  analysis  of  the  model  and 
the  actual  system. 
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IV  Example  Experiment 


Introduction 

In  order  to  demonstrate  the  accomplishment  of  two  of 
the  thesis  objectives  and  to  further  the  validation  effort, 
an  example  experiment  is  designed.  The  experiment  also  illus¬ 
trates  one  prospective  use  of  the  model  and  partially  demon¬ 
strates  the  flexibility  offered  by  SLAM  and  the  model  design 
in  analyzing  the  sensitivity  of  the  response  variable  to 
changes  in  the  control  and  stochastic  variables.  The  question 
to  be  answered  by  the  experiment  is  from  the  point  of  view 
of  the  attackers.  However,  this  is  not  the  only  point  of 
view  for  which  the  model  can  be  exercised. 

Due  to  the  highly  sensitive  and  classified  nature  of 
nuclear  targeting  parameters,  proposed  BMD  system  character¬ 
istics,  and  existing  and  proposed  RV  characteristics,  the  full 
flexibility  of  the  model  cannot  be  demonstrated  in  a  thesis 
effort.  The  values  of  all  variables  and  parameters  in  the 
example  experiment  are  based  on  values  obtained  from  unclas¬ 
sified  papers,  periodicals,  studies,  and  conversations  with 
experts  in  tne  field  (Ref  1:43-45;  2:17;  4;  5;  8;  13:12-18; 

14;  22;  23:7;  25:368;  26). 

Experimental  Design 

The  example  experiment  is  designed  to  answer  the  fol¬ 
lowing  question: 

Assuming  the  presence  of  an  improved  BMD  system  com¬ 
posed  of  an  upper  layer,  endo-atmospheric  and  lower 
layer,  point  defense  systems,  what  is  the  best  RV 


load  mix  for  a  delivery  vehicle  to  achieve  the 
greatest  probability  of  target  destruction?  The 
load  will  consist  of  new,  maneuvering  RVs  and  cur¬ 
rent  RVs. 

An  intuitive  screening  of  the  control,  stochastic, 
and  response  variables  in  Table  III,  and  preliminary  runs  of 
the  model  during  the  verification  phase,  indicate  that  three 
variables  should  be  regarded  as  being  the  most  significant 
contributing  factors  in  answering  the  above  question.  That 
is,  the  value  of  the  response  variable  is  most  sensitive  to 
the  values  of  those  three  variables.  These  three  variables 
are  RV  mix,  CEP  of  the  new  RV,  and  the  probability  of  the 
point  defense  discriminating  a  new  RV.  In  a  statistical 
experiment,  these  variables  are  referred  to  as  factors.  The 
response  variable  of  interest  is  target  area  destruction. 

The  resulting  experimental  design  chosen  for  this  experiment 
Is  a  three-factor  design,  structured  as  follows: 


ijkl 


=  y  +  +  Bj  +  4  (a8)ij  +  (ay) 


ik 


+  (SY)jk  +  (°SY)i;jk  +  eijkl  (11) 


where 


y  =  mean 

ci^,  Bj,  Yk  3  main  effects 

(08).  .,  (ay).v,  (By)ji.  *  two-factor  interaction 

1K  effects 

(aBy).  -  =  three-factor  interaction 

1JK  effect 


eijkl 


a  error  term 
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In  such  a  factorial  experiment,  all  levels  of  a  given  factor 
are  combined  with  all  levels  of  every  other  factor  (Ref  15: 
428;  21:164). 

The  levels  of  each  factor  are  then  chosen.  Three 
levels  of  each  factor  are  chosen  to  keep  the  number  of  lev¬ 
els  and  the  resulting  costs  to  a  minimum,  and  to  gain  cer¬ 
tain  analytical  advantages  (Ref  21:154).  The  three  RV  mixes 
chosen  are  three  new  RVs  and  seven  current  RVs,  five  new  RVs 
and  five  current  RVs,  and  seven  new  RVs  and  three  current 
RVs.  To  date,  the  exact  values  of  the  CEP  of  the  new  RV  and 
the  probability  of  the  point  defense  discriminating  a  new  RV 
are  unknown  and  merely  conjecture.  The  three  CEPs  chosen 
for  the  new  RV  are  200,  400,  and  600  feet.  A  CEP  of  200  feet 
Is  a  goal  in  the  design  of  the  RV.  A  CEP  of  600  feet  Is  cho¬ 
sen  since  another  design  goal  of  the  new  RV  is  that  the  accu¬ 
racy  be  at  least  as  good  as  that  of  the  current  RVs.  The 
CEP  of  400  feet  Is  chosen  as  a  mean  value.  Three  figures 
are  chosen  as  the  probability  of  the  point  defense  discrimi¬ 
nating  a  new  RV.  These  are  .5,  .3,  and  .1.  Preliminary  runs 
of  the  model  indicate  that  if  the  discrimination  probability 
of  the  point  defense  is  higher  than  .5,  no  difference  exists 
between  using  new  and  current  RVs.  The  ti ree  factors  and 
three  levels  of  each  factor  yield  the  functional  form  of  the 
experiment  depicted  in  Figure  11. 

Such  an  experiment  is  symmetrical,  all  of  the  factors 
are  set  at  the  same  number  of  levels.  A  complete  experiment 
is  desired.  That  is,  one  response  measurement  is  desired  in 


Fig  11.  Functional  Form  of  the  Experiment 

each  cell  in  Figure  11.  This  experimental  structure  is 
described  by 

N  »  qk  (12) 

O 

where 

Ng  =»  number  of  cells  in  the  experiment 
q  *  number  of  levels  of  each  factor 
k  =  number  of  factors  in  the  experiment 
In  order  to  assess  three-way  interactions,  two  response  mea¬ 
surements  are  desired  in  each  cell.  Eq  (12)  becomes 

N  =  pqk  (13) 

where 

N  =  total  number  of  computer  runs  required 
p  =  number  of  response  measurements  in  each  cell 
The  definitions  of  q  and  k  are  as  before.  Substituting  the 
appropriate  values  into  Eq  (13)  yields 
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54  =  2  x  33 


(14) 


Thus,  54  simulation  runs  are  required  (Ref  21:155-156). 

The  data  from  the  54  simulation  runs  will  be  analyzed 
using  the  analysis  of  variance  (ANOVA)  routines  from  the 
Statistical  Package  for  the  Social  Sciences  (SPSS).  The 
ANOVA  routines  will  be  used  to  determine  if  the  damage  means 
produced  by  the  RV  mixes  differ  significantly.  If  these 
means  differ  significantly,  a  Duncan's  multiple-range  test 
will  be  accomplished  using  SPSS.  This  test  divides  the  three 
RV  mix  means  into  subgroups  such  that  any  two  means  in  a  sub¬ 
group  do  not  differ  significantly.  This  test  should  indicate 
which  RV  mix  produces  a  significant  improvement  in  the  aver¬ 
age  probability  of  target  area  destruction,  if  any  (Ref  15: 
3821 . 

The  relative  differences  among  the  target  area  dam¬ 
ages  produced  by  each  mix  is  of  interest,  as  opposed  to  the 
absolute  value  of  the  target  area  damage  produced  by  each 
mix.  Kleijnen  suggests  that  correlated  sampling  is  the 
appropriate  variance  reduction  technique  to  be  used  in  this 
case.  Correlated  sampling  involves  a  comparison  of  the  mixes 
under  the  same  conditions.  In  ether  words,  the  average  tar¬ 
get  area  damage  produced  by  each  mix  is  simulated  using  the 
same  sequence  of  random  numbers.  Common  random  numbers  mean 
that  the  average  target  area  damages  are  statistically  corre¬ 
lated.  This  technique  implies  that  each  stochastic  input 
variable  has  its  own  sequence  of  random  numbers.  As  each  mix 
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is  simulated,  input  values  are  generated  using  the  same  ini¬ 
tial  values  for  the  random  number  generators  as  in  the  pre¬ 
vious  mix.  In  addition,  further  synchronization  is  required 
to  ensure  that  identical  random  numbers  are  used  for  each  RV 
in  each  mix  at  the  same  decision  point  during  the  simulation 
run.  This  variance  reduction  technique  is  implemented  to  re¬ 
duce  the  random  error  and  increase  the  precision  of  the  exper¬ 
iment  (Ref  11:200-202;  21:203-204). 

In  order  to  implement  this  technique  and  ensure  the 
synchronization  of  random  numbers  from  run  to  run,  separate 
random  number  streams  were  used  for  RV  drop-off  rate,  the 
time  for  the  RV  to  fly  through  the  endo-atmospheric  defense, 
the  time  for  the  RV  to  fly  through  the  point  defense,  and 
the  disposition  of  the  RV  in  the  endo-atmospheric,  point, 
and  target  regions.  The  SLAM  code  accomplishes  the  follow¬ 
ing: 

u. 

1.  RV  #n  gets  the  n  1  random  number  in  streams  1, 

2,  and  3,  and  gets  the  n  seven  random  numbers  in  stream  4 
for  the  first  observation  in  each  cell. 

2.  RV  #n  gets  the  (nth  +  10)fch  random  number  in 

■i»y.  a.  y. 

streams  1,  2,  and  3,  and  gets  the  (n  +  10)  seven  random 
numbers  in  stream  4  for  the  second  observation  in  each  cell. 
The  streams  are  then  reinitialized  to  generate  the  next  two 
observations  in  the  next  cell.  Random  number  stream  1  gen¬ 
erates  arriving  RVs.  Stream  2  determines  the  time  to  fly 
through  the  endo-atmospheric  zone.  Stream  3  determines  the 
time  to  fly  through  the  point  zone.  The  first  number  in  the 
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seven  from  stream  ^  determines  the  miss  distance  of  the  endo- 
atmospheric  interceptor,  the  second  determines  the  discrimi¬ 
nation  probability  of  the  endo-atmospheric  interceptor,  the 
third  determines  whether  or  not  the  endo-atmospheric  inter¬ 
ceptor  destroys  the  RV,  the  fourth  determines  the  miss  dis¬ 
tance  of  the  point  interceptor,  the  fifth  determines  the  dis¬ 
crimination  probability  of  the  point  interceptor,  the  sixth 
determines  whether  or  not  the  point  interceptor  kills  the 
RV,  and  the  seventh  determines  the  miss  distance  of  the  RV. 

The  attackers  are  assumed  to  be  ten  RVs  contained  on 
a  single  ICBM.  The  assumed  strategy  of  the  attackers  is  to 
send  the  new  RVs  through  first  in  order  to  draw  as  many  ABMs 
as  possible  and  still  retain  a  good  defense  penetration 
probability.  The  RV  drop-off  rate  is  assumed  to  be  normally 
distributed  with  a  mean  of  .5  and  a  standard  deviation  of 
.12  seconds.  The  assumed  values  of  the  new  and  current  RV 
parameters  are  as  indicated  in  Table  IV. 

The  assumed  capabilities  of  the  defenders  are  de¬ 
picted  in  Table  V.  Both  defenses  are  activated.  The  temper¬ 
ature  of  the  X-rays  from  the  endo-interceptor ' s  weapon  is 
assumed  to  be  eight  thousand  electron  volts.  In  addition, 
the  ratio  of  the  number  of  interceptors  to  the  number  of  RVs 
is  assumed  to  be  1.5  to  1. 

The  assumed  target  area  consists  of  ten  hard  targets, 
each  protected  by  one  point  interceptor.  The  sure  safe  and 
sure  kill  vulnerability  levels  for  each  target  are  assumed 
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to  be  2000  pounds  per  square  inch  (psi)  and  6000  psi,  respec¬ 
tively. 

Table  VI  reflects  the  data  generated  by  the  experi¬ 
ment.  The  data  was  obtained  by  accomplishing  three  sets  of 
simulation  runs,  obtaining  18  observations  from  each  set. 

The  SLAM  code  in  Appendix  C  was  used  to  obtain  the  18  obser¬ 
vations  in  the  first  row  of  Table  VI.  The  SLAM  code  in  Table 
VII  supplements  the  code  in  Appendix  C  for  this  experiment. 

It  should  be  noted  that  the  code  in  Appendix  C  has  the  num¬ 
ber  of  new  RVs  set  to  three  and  the  point  discrimination 
probability  set  to  .5.  To  obtain  the  data  in  rows  two  and 
three  of  Table  VI,  the  CEP  of  the  new  RV  [XX(22)]  in  the 
basic  code  is  changed  first  to  400  and  then  to  600. 

Experiment  Results  and  Conclusions 

The  results  of  a  three-way  ANOVA  are  presented  in 
Table  VIII.  These  results  indicate  that,  at  the  95%  confi¬ 
dence  level,  CEP  is  not  a  significant  main  effect.  That  is 
to  say  that  the  target  area  damage  is  not  sensitive  to  the 
new  RV  CEP.  At  this  same  level,  none  of  the  two-way  or 
three-way  interactions  are  significant.  The  ANOVA  also  indi¬ 
cates  that  there  is  a  significant  difference  in  the  average 
probability  of  target  area  destruction  produced  by  each  mix. 

Based  on  the  lack  of  significance  of  the  three-way 
interactions,  the  three-way  ANOVA  was  again  accomplished, 
forcing  the  three-way  interactions  into  the  error  term.  This 
was  accomplished  to  determine  if  the  two-way  interactions 
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table  vii 


SLAM  Supplementary  Code 


SEEBS.»(l)/N9»3(2J;HO»0(3}/NO.3(4)/NO! 

simulate; 

simulate; 

SEEDS.  IU!/N0»3(2)  /MO.  3(3)  /N0»tf(41/N0! 

iNTLC.XXUU=.3! 

simulate; 

simulate; 

SEEDSt!U!/NO.8(2)/NOi0(3)/Nfl.8(4)/NO! 

IHTLCiH(ll)».l! 

SIMULATE! 

SIMULATE! 

INTLCiH(l)*5iUm>*.5t 

SEEDS,! U)/N0»«{2)/N0.3<3}/N0.fl(4) /NO! 

simulate; 

simulate; 

SEEDS»IU»/N0»«{2)/N0,g<3)/N0»«{4)/N0! 

INTLC.XXUU-.3! 

simulate; 

simulate; 

SEEDS. fill /NO>0<2) /NOf *(3) /HOt 0(41 /*»! 
INTUSKUH*.!! 

SIBILATE! 

simulate; 

mC.n(l):7,XMIK5! 

SEEDS.  !(1) /NO. 3(2)  /NO,  !(3)/N0,  8(41  /MO! 

simulate; 

simulate; 

SEEDS.I(U/N0,!(2)/N0,fH3)/N0,!(4)/N0! 

INTLC,XX(U)=.3! 

SIMULATE! 

SIMULATE! 

$EEDS.I(1)/N0,!(2)/N0,!{3J/N0»!(4)/NQ! 

l!JTLC,XX<ll)=.H 

SIMULATE!  * 

FIN! 


INITIALIZE  SEEDS.  2  RUMS 

REINITIALIZE  SEEDS.  SET  DISCRIMINATION 
PROBABILITY  TO  .3.  2  RUNS 


REINITIALIZE  SEEDS .SET  DISCRIMINATION 
PROBABILITY  TO  .1. 


SET  NUMBER  OF  ARVS  TO  5.  SET 
DISCRIMINATION  PROBABILITY  TO  .5. 
REINITIALIZE  SEEDS.  2  RUNS 

REINITIALIZE  SEEDS.  SET  DISCRIMINATION 
PROBABILITY  TO  .3.  2  RUNS 


REINITIALIZE  SEEDS.  SET  DISCRIMINATION 
PROBABILITY  TO  .1.  2  RUNS 


SET  NUMBER  OF  ARVS  TO  7,  SET 
DISCRIMINATION  PROBABILITY  TO  .5. 
REINITIALIZE  SEEDS.  2  RUNS 

REINITIALIZE  SEEDS.  SET  DISCRIMINATION 
PROBABILITY  TO  .3.  2  RUNS 


REINITIALIZE  SEEDS.  SET  DISCRIMINATION 
PROBABILITY  TO  .1,  2  RUNS 
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changed  in  significance  and  identify  the  most  significant 
two-way  interaction.  The  results  are  depicted  in  Table  IX. 
The  degrees  of  freedom  in  the  error  term  increased  and  the 
mean  square  error  of  the  error  term  decreased.  This  indi¬ 
cates  a  desirable  result.  The  CEP  is  still  not  significant 
at  the  95%  confidence  level.  The  two-way  interactions  remain 
insignificant.  Again,  a  significant  difference  is  indicated 
in  the  target  area  destruction  produced  by  each  mix. 


TABLE  IX 

Three-Way  ANOVA  Results  -  Higher  Interactions  Confounded 


tiiitiii  ANALYSIS  OF  VARIANCE  ******** 
PROB 
BT  CEP 
MIX 
DP 


SUM  OF 

MEAN 

SIGNIF 

SOURCE  OF  VARIATION 

SQUARES 

DF 

SQUARE 

F  , 

OF  F 

MAIN  EFFECTS 

.138 

8 

.823 

8.825 

.881 

CEP 

.014 

2 

.887 

1.845 

•  .173 

MIX 

.041 

2 

.821 

5.373 

.889 

DP 

.183 

2 

.842 

18.857 

.881 

2-MAT  INTERACTIONS 

.885 

12 

.885 

1.418 

.288 

CEP  MIT 

.883 

4 

.881 

.227 

.922 

CEP  DP 

.828 

4 

.887 

1.843 

.143 

MIX  DP 

.833 

4 

.888 

2.159 

.894 

EXPLAINED 

,283 

18 

.811 

2.948 

.883 

RESIDUAL 

.134 

35 

.884 

TOTAL 

.337 

53 

.888 
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Based  on  the  resulting  higher  significance  of  the 
MIX-DP  interaction,  a  two-way  ANOVA  was  accomplished.  The 
results  are  depicted  in  Table  X.  The  degrees  of  freedom  in 
the  error  term  again  increased  and  the  mean  square  error  of 
the  error  term  remained  the  same.  This,  again,  indicates  a 
desirable  result.  The  t  way  interactions  are  still  not 
significant  at  the  95%  level.  The  target  area  damages  pro¬ 
duced  by  each  mix  remain  significantly  different.  In  all 
cases,  the  discrimination  probability  of  the  point  defense 
is  the  most  significant  frctor  in  the  probability  of  target 
area  destruction. 


TABLE  X 

Two-Way  ANOVA  Results 

hhhh  AHALYSiS  OF  VARIANCE******** 
PROB 
BT  MIX 
OP 


*********** 

************ 

*  * 

*  *  *  *  » 

*  *  *  * 

*  *  * 

SUN  OF 

MEAN 

SIGNIF 

SOURCE  OF  VARIATION 

SQUARES 

DF 

SQUARE 

F 

OF  F 

RAIN  EFFECTS 

.124 

4 

.031 

7.775 

.001 

RIX 

.041 

2 

.021 

5.147 

.010 

OP 

.083 

2 

.042 

10.402 

.001 

2-OAT  INTERACTIONS 

.033 

4 

.008 

2.069 

.101 

RIX  DP 

.033 

4 

.008 

2.069 

.101 

EXPLAINED 

.157 

8 

.020 

4.922 

.001 

RESIDUAL 

.180 

45 

.004 

TOTAL 

.337 

53 

.006 
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In  order  to  determine  and  confirm  the  recommended 


mix  to  answer  the  initial  question,  a  Duncan’s  multiple-range 


test  was  accomplished  at  the  95%  confidence  level.  The  re¬ 


sults  are  depicted  in  Table  XI.  The  results  indicate  that 


there  is  no  significant  difference  between  the  mixes  of  three 


new  RVs  and  seven  current  RVs,  and  five  new  RVs  and  five  cur¬ 


rent  RVs.  The  mix  of  seven  new  RVs  and  three  current  RVs  is 


significantly  different  from  the  other  two.  Inerefore,  based 


on  this  experiment  and  the  analysis,  a  mix  of  seven  new  RVs 


and  three  current  RVs  is  recommended. 


TABLE  XI 


Duncan's  Multiple  Range  Test  Results 


DUNCAN  PROCEDURE 
RANGES  FOR  THE  .851  LEVEL 


2.84  2.99 


THE  RANGES  ABOVE  ARE  TABULAR  VALUES. 

THE  VALUE  ACTUALLT  COMPARED  WITH  MEANIJ)-HEAN(I)  IS.. 
.9538  *  RANGE  *  SQRTU/N(I>  +  l/N(JU 


HOMOGENEOUS  SUBSETS  (SUBSETS  OF  GROUPSi  WHOSE  HIGHEST  AND  LOWEST  MEANS  DO 
NOT  DIFFER  81  MORE  THAN  THE  SHORTEST  SIGNFICANT  RANGE  FOR  A 
SUBSET  OF  THAT  SIZE) 


SUBSET  1 


GROUP  GRP  I 
MEAN  .8308 


GRP  2 
,0425 


SUBSET  2 


GRP  3 
.8943 


V  Conclusions  and  Recommendations 


Although  the  primary  goal  of  this  thesis  was  the  con¬ 
struction  of  a  simulation  model,  the  development  and  subse¬ 
quent  exercise  of  the  model  have  lead  to  some  conclusions  and 
recommendations  about  the  general  nature  of  the  RV/ABM  system. 
The  conclusions  are  based  both  on  the  objectives  of  the  thesis 
and  on  the  results  of  the  analysis  of  the  system. 

Conclusions 

The  three  objectives  of  this  thesis  have  been  satis¬ 
fied.  First,  a  verified  and  validated  simulation  model  has 
been  developed  that  will  improve  the  analysis  of  the  RV/ABM 
engagement  system.  Second,  a  "User's  Manual"  for  this  model 
has  been  provided,  by  way  of  Appendix  C,  that  instructs  the 
user  in  model  operation.  Finally,  an  illustrative  example 
is  included  that  demonstrates  both  the  flexibility  of  the 
model  and  the  complexities  of  the  actual  system.  Tne  experi¬ 
mental  design  that  is  contained  within  the  example  also  leads 
to  some  conclusions  regarding  the  statistical  significance 
of  the  factors  in  the  experiment. 

ICBM  Load  Mix.  The  first  statistically  significant 
factor  considered  is  the  load  mix  of  the  ICBM.  Assuming  the 
weapon  system  parameters  contained  in  the  model,  the  MRV  is 
more  effective  than  the  BRV  using  a  criterion  of  target  level 
damage.  This  appears  to  be  true  for  two  primary  reasons. 

The  first  is  that  the  MRV  is  more  successful  at  penetrating 


the  ballistic  missile  defenses  than  the  BRV.  The  second,  and 
less  significant,  reason  is  that  the  accuracies  assumed  for 
the  MRV  are  more  precise  than  those  assumed  for  trie  BRV, 

The  higher  penetration  success  stems  from  two  factors 
the  vulnerability  levels  of  the  RVs,  and  the  probabilities 
of  discrimination  assessed  for  the  defense.  With  increasing 
emphasis  on  achieving  higher  levels  of  nuclear  hardness,  the 
MRV,  representing  newer  technology,  would  logically  be  more 
resistant  to  nuclear  effects.  This  is  much  less  important, 
however,  than  the  lower  probability  of  discrimination  used 
for  the  MRV.  As  described  earlier,  the  probability  of  dis¬ 
crimination  is  the  summation  of  a  great  number  of  factors 
such  as  launch  and  weapon  reliabilities,  correct  object  iden¬ 
tification  and  trajectory  prediction,  and  the  probability 
that  all  defensive  systems  function  correctly  during  the 
engagement.  Most  of  these  problems  for  the  defense  are  mag¬ 
nified  when  facing  a  maneuvering  RV.  Therefore,  the  probabiL 
ity  of  discrimination  should  be  lower.  This  results  in  in¬ 
creased  success  at  defense  penetration  and  a  greater  likeli¬ 
hood  of  causing  target  damage. 

Once  at  the  target  area,  the  increased  accuracy  of 
the  MRV  appears  to  at  least  partially  offset  the  larger  yield 
of  the  BRV.  This  was  not  statistically  tested,  however. 

What  was  tested,  and  statistically  proven,  is  that  a  mix  of 
seven  MRVs  with  three  BRVs  results  in  a  statistically  signif¬ 
icant,  higher  level  of  target  damage  than  the  other  loads 
tested  (fewer  MRVs/more  BRVs). 
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MRV  Accuracy.  The  second  factor  considered  is  the 
accuracy  (CEP)  of  the  MRV.  Previous  analyses  have  identified 
the  fact  that  accuracy  is  more  important  than  yield  in  attack¬ 
ing  hard  targets  (Ref  2).  With  targets  such  as  the  hypothet¬ 
ical  ones  of  this  thesis,  hardened  to  resist  several  thousands 
of  pounds  of  overpressure  per  square  inch,  extremely  precise 
accuracy  is  required  for  a  successful  attack. 

Probability  of  Discrimination.  Perhaps  the  most  sig¬ 
nificant  result  of  this  thesis  is  the  identification  of  rhe 
criticality  of  the  probability  of  discrimination  for  the 
defense.  The  trial  runs,  mentioned  in  Chapter  III,  that  were 
used  to  test  the  assumptions  of  certain  probability  distribu¬ 
tion  functions,  revealed  a  very  important  finding.  If  the 
RV  was  successfully  identified,  tracked,  and  engaged,  it  was 
destroyed.  The  weapons  used  by  the  interceptors  are  large 
enough  that,  when  paired  with  the  precise  accuracies  of  the 
interceptors,  the  RVs  cannot  survive  the  vast  majority  of 
the  detonations.  The  significance  of  this  figure,  then, 
leads  to  the  conclusion  that  the  technology  to  permit  success¬ 
ful  penetration  to  the  target  is  as  important  as  the  tech¬ 
nology  that  permits  increased  RV  accuracy.  The  trial  runs 
that  resulted  in  forty  RV  kills  out  of  forty  engagements 
with  the  engagements  taking  place  at  altitudes  least  favor¬ 
able  to  the  defense  attest  to  this.  The  statistical  analysis 
that  verified  that  the  probability  of  discrimination  was  a 
statistically  significant  factor  in  determining  target  damage 


levels  also  is  evidence  for  this  point.  This  results,  how¬ 
ever,  in  being  able  to  go  beyond  the  maneuvering  reentry 
vehicle  to  the  conclusion  that  any  technology  that  lowers 
the  probability  of  discrimination  should  result  in  higher 
target  damage  levels.  Better  decoys  and  chaff,  electronic 
counter  measures,  and  defe  lse  suppression  are  all  options 
that  could  decrease  the  probability  of  discrimination  and 
should,  therefore,  be  investigated. 

The  implications  of  this  conclusion  are  far-reaching 
The  first  concerns  the  credibility  of  the  ICBM  deterrent. 

The  primary  concern  voiced  about  the  prelaunch  survivability 
of  the  ICBM  and  the  impact  that  it  has  on  the  credibility 
only  addresses  one  aspect  of  the  problem.  In  order  to  be  a 
credible  deterrent,  the  ICBM  and  its  RVs  must  be  capable  of 
successfully  penetrating  enemy  defenses  and  striking  their 
assigned  targets.  Without  survivability  in  both  modes,  pre- 
launch  and  penetration,  the  credibility  of  the  deterrent  is 
weakened,  and,  thus,  a  destabilizing  factor  in  the  political 
and  strategic  arena  is  introduced.  The  technology  is  at 
hand  for  the  Soviet  deployment  of  an  effective  ballistic 
missile  defense  system.  Those  who  persist  in  using  the  meta¬ 
phor  for  that  defense  of  "trying  to  hit  a  bullet  with  a  bul¬ 
let"  have  been  passed  by  the  times.  The  Soviets  clearly 
believe  the  technology  is  available  and  their  perceptions 
are  the  only  ones  that  matter.  If  they  feel  that  their 
defenses  would  be  effective  in  preventing  or  limiting  damage 
from  a  full-scale  U.3.  ICBM  strike,  then  the  ICBM  deterrent 


is  not  credible  and  the  billions  of  dollars  that  will  be  spent 
for  the  MX  missile  system  will  be  wasted. 


Recommendations 

The  intent  of  this  section  is  to  offer  recommenda¬ 
tions  for  further  study  in  the  topic  area.  It  is  certainly 
a  subject  deserving  of  more  attention  because  of  the  severe 
impact  a  credible  BMD  has  on  the  ICBM  deterrence  capability. 
The  U.S.  government  is  spending  billions  of  dollars  to  make 
the  Soviets  believe  the  MX  can  survive  an  enemy  first  strike. 
That  money  will  be  nearly  worthless  if  the  Soviets  do  not 
believe  the  RVs  the  MX  will  deliver  are  survivable  during 
penetration. 

The  first  recommendation,  therefore,  is  a  classified 
study  using  the  actual  parameters  and  capabilities  of  cur¬ 
rent  and  future  systems.  Although  this  unclassified  effort 
revealed  much  of  value,  the  actual  data  were  classified  and 
could  not  be  used.  As  a  result,  the  full  capabilities  of 
the  model  could  not  be  demonstrated.  Along  with  the  capabil¬ 
ities  of  U.S.  systems,  the  refinement  of  actual  probability 
of  discrimination  values  should  be  possible.  This  critical 
factor  is  vital  in  analyzing  BMD  impact  and  effectiveness. 

The  second  recommendation,  and  perhaps  most  appro¬ 
priate  for  an  AFIT  thesis  effort,  is  a  study  of  the  impact 
of  tactics  for  the  maneuvering  RV.  For  example,  defense 
suppression  could  have  a  very  real  impact  on  the  number  of 
subsequent  RVs  that  reach  the  target.  Suitable  targets 
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exist,  such  as  a  non-mobile  acquisition  radar.  The  effective¬ 
ness  of  such  tactics  would  certainly  be  of  interest. 

A  third,  and  related,  recommendation  is  the  evalua¬ 
tion  of  defense  tactics.  Commonly  accepted  tactics  used  to 
defeat  maneuvering  threats  such  as  netting,  or  sending  multi¬ 
ple  interceptors  after  a  single  threat,  may  not  be  required 
in  light  of  new,  highly  maneuverable  interceptors.  This 
would  reduce  the  number  of  interceptors  required  per  engage¬ 
ment,  thus  negating  an  accepted  offensive  tactic,  exhaustion 
of  defenses,  or  forcing  the  defense  to  us<~  all  interceptors 
early  in  the  battle. 

A  final  recommendation  is  the  addition  of  cost  con¬ 
siderations  in  tne  evaluation  of  the  two  competing  reentry 
vehicle  designs.  This  thesis  used  a  single  criterion  for 
measuring  effectiveness,  the  level  of  target  damage.  The 
addition  of  a  second  criterion  would  not  add  significant 
difficulty,  but  could  add  significant  insight. 
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APPENDIX  A 


X-RAY  AND  NEUTRON  FLUENCE  DERIVATIONS 


X-ray  Calculations  (Ref  4) 

Two  definitions  should  be  noted  prior  to  a  discussion 

of  the  derivations  behind  the  calculations.  The  first  of  the 

terms  to  be  defined  is  macroscopic  cross  section.  \i,  .  Macro- 

t 

scopic  cross  section  is  the  probability  of  a  reaction  by  a 
single  particle  divided  by  the  length  of  path  traveled  by  the 
particle.  Stated  differently,  it  is  the  target  atom  density 
times  the  microscopic  cross  section.  Microscopic  cross  sec¬ 
tion  is  the  target  area  of  the  particle  presented  for  inter¬ 
action  and  is  the  second  term  to  be  defined. 

The  Attenuation  Law  is  applicable  to  all  streaming 
particles.  The  law  states  that  the  decrease  in  particle 
population  over  distance  is  given  by 

dn  *  ~n(r)utdr  (15) 

where  n(r)  is  the  particle  population  at  a  distance,  r.  The 
solution  to  this  differential  equation  is 

n(r)  =  n(0)e  t  (16) 

where  n(0)  is  the  initial  particle  population.  Including 
spherical  divergence  yields 

n(r)  -  ^-0)e  —  (17) 

4ur2 

Applying  the  above  to  X-rays  as  packets  of  energy  from  a 
polyenergetic  source  yields 


F(r) 


/- 


S(hv)e  t 
^ur2 


d(hv) 


(18) 


where 

F(r)  =  X-ray  energy  fluence  at  distance  r  in  calories 
per  square  centimeter 

S(hv)  =  Planckian  source  spectrum  in  calories  per 
thousand  electron  volts 
=  macroscopic  cross  section  of  air,  which  is 
also  a  function  of  X-ray  energy  (hv) 

A  build  up  factor  (BUF)  is  included  in  the  above 
equation  due  to  the  scattering  of  the  X-rays  by  the  air  yield¬ 
ing 


F(r) 


/ 


(BUF)S(hv)e  t 
4irr* 


d(hv) 


(19) 


The  nuclear  bomb  is  a  black  body  radiator  and,  there¬ 
fore,  the  Planckian  spectrum,  represented  by  Figure  12,  des¬ 
cribes  the  source.  Planck’s  spectrum  is  defined  by  the 
function 


S(hv) 


15  (hv)3 

^TT-  e«w/lcT)  .  j 


yielding  the  density  function  in  Figure  12, 
where 


(20) 


* 


m 

M 


hv  »  energy  of  the  X-rays 

kT  =  kinetic  temperature  of  the  X-rays 


CUMULATIVE  DISTRIBUTION  FUNCTION 


Integrating  the  above  equation  yields  the  cumulative  distri¬ 
bution  function  in  Figure  12  as  follows 


G(hV)  *  1^7"  r~ ChWkT)  .  ,  «hv)  (21) 

J  o 

G(»)  =  j  "Chv7 i§y)3-~  d(hv/kT)  =  1  (22) 

Letting  x  »  hv/kT  to  normalize  the  X-ray  energy  yields 

G(x)  =  [  —Tf— -  d(x')  (23) 

*  Jo  e  -  1 

Setting  G(x)  equal  to  .1,  .2,  .3,  .4,  .5,  .6,  .7,  .8,  .9, 
and  1,  and  solving 


G  (x ) 


d(x* ) 


(24) 


for  xi  yields  xi  equal  to  1.54,  2.1,  2.58,  3.04,  3.5,  4.02, 

4.62,  5.38,  6.56,  and  10,  respectively.  This  yields  10  equal 

energy  groups  with  10$  of  the  energy  in  each  group  as  follows 

group  1  -  0-1.54 

group  2  -  1.54-2.1 

gi’Oup  3  -  2.1-2.58 

group  4  -  2.58-3.04 

group  5  -  3.04-3.5 
group  6  -  3.5-4.02 

group  7  -  4.02-4,62 

group  8  -  4.62-5.38 

group  9  -  5.38-6.56 

group  10  -  6.56-10 

Finding  the  midpoint  of  each  group  yields 
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group 

1 

- 

.77 

group 

2 

- 

1.82 

group 

3 

- 

2.34 

group 

4 

- 

2.81 

group 

5 

- 

3.27 

group 

6 

- 

3.76 

group 

7 

- 

4.32 

group 

8 

- 

5.0 

group 

9 

- 

5.97 

group 

10 

- 

8.28 

Therefore,  knowing  the  kinetic  temperature  of  the  source 
X-rays,  the  source  can  be  related  to  the  normalized  midpoints, 
x^,  as  follows 


(hv)±  =*  xt(kT)  (25) 

Due  to  the  leakage,  the  upper  limit  on  the  X-ray 
yield  is  approximately  80?  of  the  total  yield  of  the  weapon. 
Thus  the  total  X-ray  energy  possible  in  calories  per  square 
centimeter  is 


S(hv)  *  .8  (yield)  (26) 

where  the  yield  is  in  kilotons  and  one  kiloton  is  equal  to 

12  10  10 

10x  calories.  Note  that  Z  S.  =  1;  therefore,  Z  S. (.8) 

i=l  1  i=l  1 

(yield)  =  S(hv).  Thus  Si  =  1.  Note  also  that  S(hv)  is  not 
a  function  of  range. 

-wtr  -(ut/p)(pr) 

In  the  above  equation,  e  is  equal  to  e 

where  yfc/p  is  the  X-ray  attenuation  coefficient  for  air  in 

square  centimeters  per  gram.  This  coefficient  is  a  function 

of  X-ray  energy,  hv,  and  is  different  for  each  of  the  10 

energy  groups  just  mentioned.  The  attenuation  coefficient 


per 


for  each  group  as  a  function  of  hv^  in  square  centimeters 
gram  is  given  by 


(ut/p)i  »  -.00135^  +  l9.756H/(hv)1  -  46l.763/(hv)| 


+  6680.02/(hv)’  -  3497.36/(hv)J 


+  907.358/(hv)5  (27) 


where 

p  =*  density  of  the  air  in  kilograms  per  cubic  meter 
r  =  range  from  burst  to  target  in  kilometers 
The  density  of  the  air  times  the  range  is  known  as 
the  mass  integral.  The  mass  integral  is  a  function  of  the 
altitude  of  the  burst,  the  altitude  of  the  target,  and  the 
range  between  the  target  and  the  burst.  If  the  burst  and 
the  target  are  at  the  same  altitude,  the  mass  Integral  is  a 
constant  for  all  groups  given  ar.  altitude  and  a  range  of  the 
engagement,  and  is  in  units  of  grams  per  square  centimeter. 

BUF  is  a  function  of  ut/p  and  is  different  for  each 
of  the  10  energy  groups.  BUF  is  unitless  and  is  given  by 


BUF1  =  Ale 


c1(wt/p)1(pr) 


+  A2e 


c2(vt/p)i(pr) 


(28) 


where  Al,  c^,  and  c2  are  dependent  on  hv^  and  are  given  in 
Tabel  XII.  In  addition,  A2  =  1  -  Al,  The  above  equation  for 
BUF  is  a  fit  for  Figure  13. 

Knowing  that 


TABLE  XII 


Constants  for  Empirical  Build-up  Factor  Equation  (Ref  4) 


Energy  in  keV 

A1 

*2 

C1 

C2 

Maximum#  difference 
below  8  HFP 

.  12 

-0.227 

1.227 

-0.400 

0.000 

0.35  * 

14 

-0.370 

1.370 

-0.400 

0.000 

0.89  % 

16 

-O.323 

1.323 

-0.680 

0.020 

1.78  * 

18 

-0.634 

1.634 

-0.460 

0.020 

0.91  % 

20 

-1.072 

2.072 

-O.36O 

0.020 

2.71  % 

22 

-1.048 

2.048 

-0.460 

0.040 

2.90  % 

24 

-1.740 

2.748 

-0.340 

0.040 

0.75  % 

26 

-2.673 

3.673 

-0.260 

0.040 

1.85  % 

28 

-2.664 

3.664 

-0.300 

0.060 

3*87  % 

30 

-6.038 

7.038 

-0.140 

0.040 

2.90  % 

32 

-8.805 

9.805 

-0.100 

0.040 

1.45  % 

34 

-8.504 

9.504 

-0.100 

0.060 

1*63# 

36 

-75.83 

76.83 

0.000 

0.020 

0.97  % 

38 

-20.03 

21.03 

-0.020 

0.060 

3-77% 

40 

-16.94 

17.94 

-0.020 

0.080 

6.93* 

45 

14.59 

-13.59 

0.120 

-0.020 

13-7* 

50 

11.31 

-IO.31 

0.160 

-0.040 

21.1  % 

55 

109.2 

-108.2 

0.120 

0.100 

18.0  % 

60 

-U.05 

12.05 

0.000 

0.200 

26*6  % 

70 

-114.1 

115.1 

0.140 

0.160 

25*7  * 

80 

-113.1 

114.1 

0.160 

0.180 

27.6  % 

too 

-10.93 

11.93 

0.060 

0.260 

36.0  % 

120 

-8.153 

9.153 

0.020 

0.280 

39.1  % 

150 

13.14 

-12.14 

0.260 

0.100 

36.3  3* 

200 

-88.92 

89.92 

0.180 

0.200 

33.1  % 

250 

-6.308 

7.308 

0.000 

0.260 

39.9  % 

300 

19.89 

-18.09 

0.200 

0.120 

33.1  * 

350 

72.37 

-71.37 

0.160 

0.140 

31-2  # 

400 

-6.O63 

7.063 

0.000 

0.220 

35.1  % 

500 

16.29 

-15.29 

0.160 

0.080 

28.8  % 

600 

-57.58 

57.56 

0.100 

0.120 

25.7  % 

750 

-17.20 

18.20 

0.060 

89 

0.120 

22.6  % 

0 « 00  3  on 

£1  stance  from 


Flg  13 ’  £nergy  Build 


LU^n-frYfV^ 

3rS  ln  H°»°eeneou a  Air  rB.,  7 


d(hv) 


(29) 


F(r) 


~vtr 

(BUF)S(hv)e  z 


4irr2 


and  transforming  F(r)  into  a  discrete  function  with 


10 

BUF  =  z  BUF,  (30) 

i=l 


10 

S(hv)  =  z  S,  ( .  8 ) (yield)  (31) 

i=l  1 


where  *  .1  and 


-y  r  -(u4./p)(pr)  10  -(y. /p ), (pr) 

et=ec  =  Z  e  z  1  (32) 

i=l 


yields 


F(r) 


.8 (yield) 
4irr2 


10  -(u,./p  ) .  (pr) 

l  ( . 1) (BUF) . (e  °  1  ) 

i=l  x 


ca]./cm2 


(33) 


This  formulation  is  good  for  X-rays  with  kinetic  temperatures 
which  are  greater  than  or  equal  to  a  few  thousand  electron 
volts  because  it  ignores  scatter. 

User  Function  6  uses  the  above  X-ray  equations.  Two 
assumptions  were  made  concerning  the  X-ray  fluence  on  the 
target.  The  first  of  these  is  that  the  atmosphere  is  homo¬ 
geneous.  The  second  is  that  X-rays  are  considered  damaging 
above  70,000  feet,  only.  Given  the  yield  in  kilotons,  the 
temperature  of  the  source  X-rays  in  kT,  and  the  range  from 
the  burst  to  target  in  feet,  the  FORTRAN  code  in  User  Function 


;s 


6  first  calculates  the  yield  in  calories  as  follows 

yield(cal)  -  yield(kilotons) (1012  calories/kiloton)  (34) 

The  range  is  then  converted  to  kilometers  as  follows 

r(km)  =  r(ft)(. 00030480  km/ft)  (35) 

The  mass  integral  is  then  computed  by 

pr(gm/cm2)  *  p (kg/m3 )r (km) ( 103  gm/kg)(10^  cm/km) ( 10“^  m3/cm3) 

(36) 

=  p (kg/m3 )r(km) ( 102 )  (37) 

where  p  is  the  density  of  the  air  at  the  intercept  altitude. 

For  each  of  the  10  groups,  hv^  is  computed  using  Eq  (25),  and 

( /p ) ^  is  computed  using  Eq  (27).  An  iterative  routine  is 

then  accomplished  to  find  the  nearest  hv  in  Table  XII  to  hv^ 

The  hv  value  obtained  is  then  used  to  obtain  Al,  A2,  c^}  and 

Cg.  BUF.^  is  then  computed  using  Eq  (28).  For  each  of  the 

-(ut/p)i(pr) 

ten  groups,  a  value  is  computed  by  (,l)(BUF)^e  1 

As  each  value  is  obtained  for  a  group,  they  are  'summed  to 
obtain  the  sum  for  all  10  groups.  The  range  is  then  con¬ 
verted  to  centimeters  as  follows 

r(cm)  =  r(km)(10^  cm/km)  (38) 


Finally,  the  total  X-ray  fluence  on  target  is  computed  using 
Eq  (33).  It  should  be  noted  that  if  hv^  is  less  than  12, 


BUFi  is  set  equal  to  1. 
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Neutron  Calculations  (Ref  4) 


In  contrast  with  X-rays,  the  Attenuation  Law  is  not 
applicable  to  neutrons.  Therefore,  Eq  (39)  represents 

p(v)  _  3. 16E23  (neutrons/klloton)yield(kilotons)S(x)  ^g) 

4Trr(cm) 2 

the  neutron  fluence  on  the  target  in  neutrons  per  square 
centimeter.  The  above  equation  assumes  a  thermonuclear  yield 
which  produces  3.16E23  neutrons  per  kiloton  of  yield. 

Absorption  of  neutrons  by  the  atmosphere  is  not  likely. 
Therefore,  scatter  dominates.  The  repeated  scattering  and 
resulting  energy  loss  implies  that  the  average  energy  of  the 
neutrons  decreases  as  time  and  distance  form  the  burst  in¬ 
crease.  The  procedure,  then,  for  determining  the  neutron 
< 

fluence  per  source  neutron,  S(x),  at  a  particular  range  is 
to  find  the  fluence  in  homogeneous  air  and  then  apply  the 
result  to  every  point  in  the  real  atmosphere  which  lies  on 
the  same  mass  range  from  the  burst  point.  This  procedure 
assumes  that  spherical  divergence  and  neutron  attenuation 
are  separable,  which  they  are  not.  Sizeable  error  results 
for  altitudes  from  5  to  25  kilometers.  However,  the  code 
which  is  currently  popular  uses  this  procedure  and  the  re¬ 
sults  are  good. 

S(x)  is,  then,  a  function  of  the  mass  integral,  pr. 

Once  again,  the  density  of  the  air,  p,  varies  with  altitude. 

If  the  burst  and  the  target  are  coaltitude,  the  mass  integral 
is  determined  given  the  altitude  and  the  range  of  the 


engagement,  and  is  in  units  of  grams  per  square  centimeter. 
u(x)  as  a  function  of  the  mass  integral  is  given  by 

S(x)  *  exp[A  +  B(pr)  +  C(pr)^  +  D(pr)^’^  +  E(pr)*3 

+  F(pr)1/3  +  Gln(pr) ]  (40) 

The  variables  A  through  G  are  determined  from  Table  XIII  for 
the  neutron  fluence  from  a  thermonuclear  source.  The  above 
equation  is  an  empirical  fit  for  Figure  14. 

User  functions  6  and  7  use  the  above  neutron  equa¬ 
tions.  An  additional  assumption  made  in  the  above  calcula¬ 
tions  Is  that  the  atmosphere  is  homogeneous  and  exponential. 
Given  the  yield  of  the  weapon  in  kilotons  and  the  range  from 
the  burst  to  the  target  in  feet,  the  FORTRAN  code  in  User 
Functions  6  and  7  first  converts  the  range  to  kilometers 
using  Eq  (35).  The  mass  integral  is  then  computed  using  Eq 
(37).  Next,  S(x)  is  computed  using 

S(x)  *  expC-6.775  +  .005296(pr)  -  . 0000053^6(pr) 2 

-  .0002l468(pr)1*5  -  3.82l4(pr)*5 
+  10.875(pr)1/3  -  1.39751n(pr)](4l) 

The  range  is  then  converted  to  centimeters  using  Eq  (38). 
Finally,  the  total  neutron  fluence  on  target  is  computed 
using  Eq  (39). 

Overpressure  Calculations  (Ref  4) 

The  miss  distance  of  the  weapon  is  scaled  down  to  a 
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yield  of  one  kiloton,  the  reference  for  graphs  and  empirical 
calculations  by 

range  =  range/yield  (42) 

The  height  of  the  burst  and  the  target  are  assumed  to  be 
zero  at  sea  level.  User  Function  8  uses  an  empirical  fit 
for  Figure  15  to  compute  the  overpressure  placed  on  the  tar¬ 
get  by  an  RV. 


400 


15.  Peak  Overpressures  on  the  Ground  for  a  1  kt  Burst 
(Ref  9:111) 


APPENDIX  B 

CONTINUOUS  DAMAGE  FUNCTIONS  AND  WEAPON  RADIUS  DERIVATIONS 


Damage  Functions  (Ref  4;  25:346-362) 

Continuous  damage  functions  were  defined  and  used  to 
evaluate  the  X-ray  fluence,  neutron  fluence,  and  overpressure 
on  the  target.  For  X-ray  and  neutron  fluences,  the  functions 
are  based  on  the  ranges  of  the  damage  mechanism.  For  over¬ 
pressure,  the  function  is  based  on  the  intensity  of  the  dam¬ 
age  mechanism.  In  all  cases,  the  damage  function  is  based 
on  a  lognormal  probability  distribution  function  where 
Pd  =  probability  of  damage 
R  =  range 

I  *  intensity 

R  ej  *  range  for  50% 

I  3  intensity  for  50?  'J?d 
The  lognormal  function  is  given  by 


f(x) 


-i  (-iSVLJL-) 


ax  y  2v 


(43) 


where  y  and  a  are  for  the  corresponding  normal  distribution. 

Considering  intensity  first,  the  cumulative  function 
is  given  by 

2 

(44) 


V1’  -f 


ol 


V27 


Letting  z  =  (lnl  -  y)/o  yields 


lnl  3  za  +  y 


(45) 
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I  =  e 


(zo  +  y) 


(46) 


When  I  a  Q,  z  =  -®,  and  dz  =  dl/al.  Therefore, 


V2> 


The  sure  kill  intensity,  I  is  defined  as  the  inten¬ 
sity  at  which  the  target  is  Killed  98/5  of  the  time.  The  sure 

safe  intensity,  I__,  is  defined  as  the  intensity  at  which  the 

s  s 

target  is  killed  2)5  of  the  time.  Thus, 


Prom  standard  normal  tables, 
z  *  -2.054  when  Pd(z)  *  *02. 

2.054  * 

-2.054  = 

or 

2.054a  = 
-2.054a  = 


=  2.054  when  Pd(z)  31  .98  and 
Therefore, 


lnIsk  -  “ 
a 

(50) 

lnIss  ~  * 
a 

(51) 

lnIsk  -  » 

(52) 

lnIss  -  “ 

(53) 

L-  @Qf 


Solving  Eqs  (52)  and  (53)  simultaneously  for  y  yields 


i 

0  =  lnIsk  +  lnIss  -  2u 

(54) 

1 

2w  -  lnlsk  +  lnlgs 

(55) 

-  ■  ‘“V..-1 

(56) 

Solving  Eqs  (52)  and  (53)  simultaneously  for  a  yields 


2(2.054)a  =  lnlsk  -  lnlss 


(58) 


Thus,  knowing  the  intensity  for  sure  kill,  the  intensity  for 
sure  safe,  and  the  intensity  on  target 


lnl  -  y 


where  u  and  a  are  computed  from  Eqs  (56)  and  (58),  respec¬ 
tively. 

Having  computed  z,  the  following  equation  can  be  used 


to  approximate  ( z ) 


Pd(z)  =  .5(1  +  .19685^z  +  .115194z2  +  .000344z3 


h  -U 

+  .019527Z  )  * 


If  z  is  less  than  0,  z  in  Eq  (60)  is  set  to  the  absolute  val¬ 
ue  of  the  computed  z.  If  z  is  greater  than  or  equal  to  0, 


Pd(z)  =  1  -  Eq  (60) 
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User  Function  8  uses  this  intensity  approach  to  cal¬ 
culate  the  probability  of  damage  of  an  RV  against  a  target 
where  I  *  2000  pounds  per  square  inch  and  I  .  ®  6000  pounds 
per  square  inch. 

A  related,  but  different  approach,  was  used  to  calcu¬ 
late  the  probability  of  damage  of  the  ABM  against  the  RV.  A 
concept  known  as  weapon  radius  (WR)  was  used.  Weapon  radius 
is  the  "cookie  iutter"  radius  that  will  destroy  the  same  num¬ 
ber  of  targets  as  does  the  actual  probability  distribution 
function.  The  weapon  radius  depends  on  the  Igs  and  I  ^  of 
the  target  as  well  as  the  yield  of  the  weapon.  The  weapon 
must  burst/impact  within  the  weapon  radius  to  kill  the  target. 
Weapon  radius  calculations  assume  homogeneously  distributed. 
Identical  targets  and  a  given  probability  distribution  of 
damage,  Pd(R). 

Range,  rather  than  intensity,  is  used  in  developing 
the  weapon  radius  where 

it  (WR) 2  *  J  Pd(R)2irRdR  (62) 

(WR)2  =  /  P  (R)2RdR  (63) 

J  o  a 

Integrating  Eq  (63)  by  parts,  letting  u  *  Pd(R)  and 
dv  =  2RdR,  yields 


d(P .  (R) ) 


du  =* — hr 

(64) 

v  =  R2 

(65) 
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Therefore . 


(**)■  "  CR2Pd(H)l  "-J^  R2  dRr 


rl 


-  d(Pd(R)) 


dR 


(66) 


However, 


0  =  [R2Pd(R)]“ 


(67) 


Therefore, 


d(P, (R) ) 

(TO)2  »  |  R2  - is -  dR 


L 


dR 


(68) 


Again  basing  the  damage  function  on  a  lognormal  prob¬ 
ability  distribution  yields 


(WR) 1 


r 


i 


_j(OnR^) 


cR  y  2ir 


dR  (69) 


Letting  z  *  (InR  -  y)/a  yields 


InR  *  za  +  u 


R  *  e(za  +  y) 


(7u) 

(71) 


When  R  =  0 ,  z  =  -  «  ,  and  dz  =  dR/Rc.  Therefore, 


(TO) 2 

(TO)2 


L 


V2zo  +  -  e-***dz 


i2“  -  e'!z‘ 


(72) 

(73) 


(WR) 2  « 


s2y 


r  to 

J- 


1  e-i(z2  -  4zo)dz 


y/27 


(74) 


(WR)2  =  e 


(WR)2  =  e 


2y 


[”1  -  4za  +  4a2)dz  (?5) 

L  v 2* 


(2y 


+  2°2>  <%*■,,.  e^<2  -  2c)2dz 

A.  V77 


(76) 


Letting  w  »  z  -  2a  and  dw  *  dz  yields 


(WR)*  =  e2(“  +  oI>  f  —i -  e-Jw2 

>ico  ,/ 2 IT 


dw  (77) 


However, 


■  ■  />■ 


(78) 


Therefore, 


(WR)2  ■  e2^w  +  ff2) 


WR  =  e(w  +  °2) 


(79) 

(80) 


Since 


Vr)  - 


/ 


R 


1  .i(-3aSrJL-)< 


crR  Vi7 


dR  (81) 


y  and  a  can  be  expressed  as 'functions  of  the  range  of 'sure 
kill,  R  .  ,  and  the  range  of  sure  safe,  R  ,  just  as  in  the 

SK  S  5 

intensity  case  previously  developed.  This  yields 


105 


*  -  •51n(IW,ss> 


(82) 


a 


1 

2(2.034) 


ln(IVW 


(83) 


Knowing  the  intensity  for  sure  kill  and  the  intensity 
for  sure  safe,  the  X-ray  and  neutron  calculation  routines  in 
User  Functions  6  and  7  are  accomplished  for  different  ranges 
to  determine  the  ranges  at  which  the  sure  safe  and  sure  kill 
intensities  are  realized.  Having  Rss  and  Rsk,  WR  can  be 
computed  for  the  particular  encounter  using  Eq  (80). 

Prior  to  further  development  of  this  damage  function, 
four  assumptions  must  be  made  explicit.  The  first  of  these 
is  that  the  target  is  in  the  burst  plane.  The  second  is  that 
the  height  of  the  burst  is  controllable  such  that  a  *  0. 

The  third  is  that  variations  in  missing  in  the  x  and  y  direc¬ 
tions  are  equal  such  that  a  =»  o  .  Finally,  the  distributions 

x  y 

in  the  x  and  y  directions  are  normal.  Therefore, 


f(x) 


f  (y)  * 


l 

5yV^* 


(84) 


(85) 


Thus,  f(x,y)  for  a  differential  area  dA  is  given  by 

if  **  +  y2  ) 

f(x,y)  *  — ~ —  e  a  (86) 

o22* 
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Transforming  into  the  circular  normal  function  yields 

Z 


i  -i(Hr) 

f(r>9)  ■  snr  e  * 


(87) 


The  cumulative  function  is  given  by 


rdrde 


(88) 


Integrating  with  respect  to  0  yields 


F(r)  = 


■/+ 


rdr 


F(r) 


(89) 


(90) 


Circular  error  probable  (CEP)  by  definition  is  the  point  at 
which  the  cumulative  function  is  equal  to  .5. 

Therefore, 


.5-j 


CEP 


-*(-§-) 


dr 


.5 


1  -  e  W  ' 


.5  =  e  V  a  ' 


(91) 


(92) 

(93) 


In. 5  =  - 


i  CEP2 
4  — tr— 


-21n. 5  =  - 


CEP! 


(9^) 

(95) 
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M  r 


CEP2  =«  (21n2)o2 


CEP  »  v 21n2  a 


21n2 


Substituting  Eq  (98)  into  Eq  (90)  yields 


F(r)  =  1  -  e 


r  V  21n2 


F(r)  =  1  -  e 


i  (  (21n2)r2  \ 
v  ceF5 - ) 


F(r)  -  1  -  e 


-ln2 


Hri 


F(r) 


=  1  -  e~*693(  CEF') 


(100) 


(101) 


(102) 


The  single  shot  probability  of  kill  (Pssk)  can  be 
computed  if  the  weapon  radius  and  CEP  of  the  weapon  are  known 
by  substituting  WR  into  Eq  (102)  for  r,  yielding 


P  *  1  -  e 
rssk  1  e 


-wUk-Y 


(103) 


In  other  words,  Pssk  is  the  probability  that  the  weapon  im¬ 
pacts  within  radius  WR  if  the  CEP  of  the  weapon  is  known. 
This  weapon  radius  rationale  is  used  in  User  Functions  6  and 


Aiming  Errors  (Ref  25:3^6-362) 


As  developed  previously,  the  aiming  errors  in  the  xy 


plane  are  assumed  to  be  Independent  and  normally  distributed 
with  equal  variance  and  zero  mean.  Therefore,  the  distribu¬ 
tion  of  miss  distances  is  given  by  the  Rayleigh  distribution 
or  circular  normal  function,  Eq  (87).  The  cumulative  func¬ 
tion  as  a  function  of  CEP  is  given  by  Eq  (102).  Solving  Eq 
(102)  for  r  yields 


CEP2ln[ 1  -  F(r) ] 
r - T7593 - 


(104) 


The  cumulative  distribution  function,  F(r),  is  uniformly  dis¬ 
tributed  between  zero  and  one.  User  Functions  6,  7,  and  8 
make  use  of  this  fact  to  compute  ABM  and  RV  miss  distances. 

On  any  given  encounter,  a  uniformly  distributed  random  num¬ 
ber  between  zero  and  one  is  drawn  and  substituted  into  Eq 
(104)  for  F(r).  Based  on  the  CEP  of  the  weapon,  the  miss 
distance,  r,  for  that  encounter  may  be  computed. 


f- 

I 


109 


Introduction 


This  appendix  contains  information  for  the  person 
using  the  RV/ABM  engagement  model  developed  in  this  thesis. 

A  knowledge  of  the  simulation  language,  SLAM,  while  helpful, 
is  not  required.  The  appendix  shows  the  job  control  lan¬ 
guage  required  for  the  CDC6600,  at  AFIT,  and  the  methods  of 
changing  the  weapon  system  parameters  for  the  offense  and 
defense.  Little  adjustment  to  the  actual  network  is  illu¬ 
strated  because  the  possibilities  are  endless  once  the  basic 
structure  is  modified.  The  design  of  the  model  is  such  that 
there  is  sufficient  flexibility  to  model  a  wide  variety  of 
situations  within  the  existing  structure.  However,  the  user 
must  have  access  to  the  classified  information  necessary  to 
take  full  advantage  of  the  capabilities  of  the  model. 

The  model  can  be  utilised  to  answer  and  analyze  ques 
tions  posed  by  both  the  offense  and  the  defense.  The  model 
is  capable  of  handling  an  analysis  of  situations  involving 
one  RV  type  or  two  RV  types,  delivered  to  the  target  in  any 
numbers  and  mixes  desired.  The  type  of  defense  can  range 
fr*om  none  at  all  to  a  layered  terminal  defense,  composed  of 
an  endo-atmospheric  and  point  defense  system.  The  target 
area  can  consist  of  any  fixed,  homogeneous  target  set. 

The  SLAM  portion  of  the  model  is  developed  in  six 
segments.  These  segments  are  discussed  in  the  following 
six  sections.  Each  section  contains  a  briei  description  of 
that  segment’s  operation,  the  SLAM  network  structure,  and 
the  SLAM  computer  code  for  that  segment  with  a  line  by  line 


interpretation  of  the  SLAi'l  computer  code.  The  SLAM  code 
Included  here  is  the  code  for  the  example  experiment  in  Chap¬ 
ter  IV.  Following  these  six  sections  is  a  listing  of  the 
commented  FORTRAN  user  functions.  The  subroutines  are  called 
and  controlled  b„-  the  SLAM  code.  For  a  more  complete  dis¬ 
cussion  of  the  user  functions,  refer  to  Appendices  A  and  B. 
Table  XIV  is  a  key  to  the  attributes  and  variables  to  be 
referred  to  in  the  following  discussion.  Any  further  changes 
beyond  those  indicated  would  require  alterations  of  the  SLAM 
network  and/or  supporting  subroutines.  With  the  aid  of  this 
appendix  and  the  example  experiment  in  Chapter  IV,  the  user 
should  be  able  to  analyze  a  wide  variety  of  RV/ABM  issues 

Arriving  RVs 

The  first  node,  START,  simulates  the  normally  distri¬ 
buted  release  of  a  specific  number  of  entities,  RVs  and 
decoys,  from  the  post-boost  vehicle.  Each  entity  is  released 
as  if  it  was  an  RV.  For  example.  If  a  decoy  weighs  10  pounds 
and  the  RV  weighs  100  pounds,  one  RV  position  could  be  used 
for  10  decoys.  Assuming  the  post-boost  vehicle  capacity  is 
10  RVs,  the  post-boost  vehicle  could  carry  10  RVs,  9  RVs  and 
10  decoys,  8  RVs  and  20  decoys,  etc.  The  number  of  entitles 
in  each  case,  however,  is  10. 

At  node  Nl,  the  entitles  are  counted.  Depending  on 


the  entity  count,  the  entity  is  then  designated  as  an  RV  of 
the  first  type,  an  ARV,  at  node  ARV,  an  RV  of  the  second  type, 
a  BRV,  at  node  BRV,  or  a  set  of  decoys  at  node  N2.  All  RVs 


TABLE  XIV 


Attribute  and  Global  Variable  Key 

ATTRIBUTE  KEY 

1  -  MARK  TIME 

2  -  RV  TYPE  (1  -ARV.  2  »  BRV) 

3  -  TIME  TO  FIT  150.088  FT  TO  70*8*0  FT  THRU  ENDO  DEFEfiSE  ZONE 

FROM  USERF  1 

4  =  TIME  TO  INTERCEPT  IN  ENBO  ZONE  FROM  USERF  2 

3  =  TINE  TO  FIT  78i888  FT  TO  18*888  FT  THRU  POINT  DEFENSE  ZONE 
FROM  USERF  3 

6  =  TINE  TO  INTERCEPT  IN  POINT  ZONE  FROM  USERF  4 

7  -  PROBABILITY  OF  ABN  KILLING  RV  IN  ENDO  OR  POINT  ZONES 

FROM  USERF  6  Ok  7 

3  *  14  =  UNIFORM  RANDAN  NO  (0  TO  1)  TO  ACHIEVE  CORRELATED  SAMPLING 

GLOBAL  VARIABLE  KET 
XX<1)  -  NO  Or  ARVS 

IX <21  -  ARV  REENTRY  ANCLE  IN  DEGREES 
XX  (3)  =  BRV  REENTRY  ANGLE  IN  DECREES 

XX (4)  -  SWITCH  fOR  ENDO  DEFENSE  (0  =  ENDO  PRESENT.  1  =  ENDO  ABSENT! 

XXloJ  -  NOT  USED  (INITIALIZED  TO  11 

XX (6)  -  PROBABILITY  THAT  ENDO  DETECTS  ARV 

XX(7)  *  VELOCITT  OF  RV  IN  FT/SEC 

XX (8)  *  PROBABILITY  THAT  ENDO  DETECTS  BRV 

XX (9)  =  TEMPERATURE  OF  X-RAYS  FROM  ENDO  ABN  IN  KEV 

XX 110)  *  SWITCH  FOR  POINT  DEFENSE  (0  =  POINT  PRESENT.  I  *  POINT  ABSENT) 

XX(ll)  *  PROBABILITY  T^AT  POINT  DETECTS  ARV 

XX (12)  -  SURE  SAFE  NEUTRON  INTENSITY  FOR  ARV  IN  NEUTR0NS/CH2 

XX (13)  *  PROBABILITY  THAT  POINT  DETECTS  BRV 

XX (14)  =  SURE  KILL  NEUTRON  INTENSITY  FOR  ARV  IN  NEUTR0NS/CN2 

XX (15)  -  PROBABILITY  OF  TARGET  DESTRUCTION  (INITIALIZED  TO  8) 

XX  (16)  :  NOT  USED  (INITIALIZED  TO  3) 

XX (17)  =  YIELD  OF  ENDO  INTERCEPTOR  IN  KILOTONS 

XX (18)  :  CEP  OF  ENDO  INTERCEPTOR  IN  FEET 

XX (19)  *  YIELD  OF  POINT  INTERCEPTOR  IN  KILOTONS 

XX (20)  =  CEP  OF  POINT  INTERCEPTOR  IN  FEET 

XX (21)  -  YIELD  OF  ARV  IN  FEET 

XX (22)  -  CEP  OF  ARV  IN  FEET 

XX (23)  =  YIELD  OF  BRV  IN  KILOTONS 

XX (24)  =  CEP  OF  BRV  IN  FEET 

XX (25)  =  COUNTER  FOR  GENERATED  RVS  (INITIALIZED  TO  0! 

I  XX (26)  =  TOTAL  NO  OF  RVS  ON  BOOSTER 
i  XX127)  :  SURE  SAFE  X-RAY  INTENSITY  FOR  ARV  IN  CAL0RIES/CM2 

!  XX (23)  *  SURE  KILL  X-RAY  INTENSITY  FOR  ARV  IN  CAL0RIES/CN2 

!  XX (29)  =  SURE  SAFE  NEUTRON  INTENSITY  FOR  BRV  IN  NEUTRGNS/CH2 

!  XX (31)  *  SURE  KILL  NEUTRON  INTENSITY  FGR  BRV  IN  NEUTR0NS/CH2 

XX (31)  =  SURE  SAFE  X-RAY  INTENSITY  FOR  BRV  IN  CAL0RIES/CN2 
XX (32)  =  SURE  KILL  X-RAY  INTENSITY  FOR  8RV  IN  CAL0RIES/CH2 
XX (33)  -  SURE  SAFE  OVERPRESSURE  FOR  TARGET 
;  XX (34)  =  SURE  KILL  OVEf PRESSURE  FOR  TARGET 
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of  the  first  type  are  introduced  into  the  system  prior  to 
release  of  any  RVs  of  the  second  type.  At  node  N2,  the  de¬ 
coys  are  assumed  to  be  stripped  away  and  do  not  enter  the 
remainder  of  the  system. 

At  node  N3,  each  RV  is  assigned  timing  information 
from  the  FORTRAN  subroutines  as  attributes.  At  node  N4,  each 
RV  i3  subsequently  assigned  a  series  of  uniform  random  num¬ 
bers  between  zero  and  one  in  order  to  achieve  correlated 
sampling  discussed  in  Chapter  IV. 

A  decision  is  made  at  node  SWE.  If  an  endo  defense 
is  not  present,  the  RV  continues  on  to  the  point  defense 
switch.  If  an  endo  defense  is  present,  the  RV  continues  to 
node  N5  where  another  decision  is  made.  If  there  are  no 
endo-interceptors  available,  the  RV  continues  on  to  the  po3nt 
defense  switch.  If  endo-interceptors  are  available,  the  RV 
is  engaged  by  the  endo  defense. 

The  line  of  code  requiring  user  change  in  this  seg¬ 
ment  is 

START  CREATE , RNORM (a,b,l),0,l,c,l; 

The  user  enters  the  mean  of  the  normally  distributed  RV  re¬ 
lease  time  for  a.  The  standard  deviation  of  these  release 
times  is  entered  for  b.  Th<-  user  enters  the  total  number  of 
entities,  RVs  and  decoys,  to  be  released  for  c. 
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Endo  Defense 


As  an  RV  reaches  node  END,  it  is  assigned  an  endo- 
interceptor  and  continues  on  to  node  N6  where  the  intercept 
time  is  reached.  At  node  N6,  the  RV  is  assigned  a  probabil¬ 
ity  of  being  destroyed  by  the  defense  from  the  FORTRAN  sub¬ 
routines  as  an  attribute.  Depending  on  RV  type,  the  RV  then 
continues  to  node  N7  or  N8.  The  network  following  these  two 
nodes  is  the  same  except  for  the  statistics  collected  at 
nodes  N13  and  Nl4,  respectively.  Following  N7,  the  number 
of  ARVs  destroyed  is  collected.  Following  node  N8,  the  num¬ 
ber  of  BRVs  destroyed  is  collected. 

Since  the  networks  are  the  same,  only  that  following 
node  N7  will  be  discussed.  At  node  N7,  a  decision  is  made 
as  to  whether  or  not  the  RV  is  discriminated.  If  discrimina¬ 
tion  does  not  occur,  the  RV  continues  on  to  the  point  defense 
switch  in  the  time  remaining  to  fly  through  the  endo  defense 
zone.  If  the  RV  Is  correctly  discriminated,  a  decision  is 
made  at  node  N9  as  to  whether  or  not  the  RV  is  destroyed. 

If  the  RV  is  destroyed,  the  number  of  destroyed  RVs  of  that 
type  is  increased  by  one.  If  the  RV  is  not  destroyed,  it 
continues  on  to  the  point  defense  switch.  All  of  the  deci¬ 
sions  are  based  on  the  outcome  of  a  random  number  draw. 
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Point  Defense  Switch 

As  an  RV  reaches  node  SV/P,  a  decision  is  made  as  to 
whether  or  not  a  point  defense  exists.  If  the  defense  does 
not  exist,  the  RV  continues  on  to  the  target.  If  a  defense 
exists,  the  RV  continues  to  node  N15  where  another  decision 
is  made.  If  there  are  no  point  interceptors  available,  the 
RV  continues  on  to  the  target.  If  point  interceptors  are 
available,  the  RV  is  engaged  by  the  point  defense. 


Point  Defense  Switch 


Point  Defense 


The  endo  defense  discussion  applies  here.  The  only 
differences  are  that  point  interceptors  are  assigned,  attri 
bute  values  differ,  and  different  FORTRAN  subroutines  are 
called.  The  network  logic  and  flow  and  the  decisions  are 
the  same. 


Point  Defense 
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As  a  surviving  RV  reaches  node  TGT,  the  RV  continues 
on  to  the  target  in  the  time  it  takes  to  fly  the  last  10,000 
feet  to  the  target.  At  node  N25,  the  average  probability  of 
target  area  destruction  is  accumulated  and  updated  as  each 
RV  reaches  its  target.  The  number  of  surviving  RVs  is  col¬ 
lected  by  type.  The  simulation,  one  run,  is  then  terminated. 


SLAM  Support  Code 

There  is  no  network  structure  for  the  SLAM  support 
code.  The  first  requirements  for  a  simulation  run  are  those 
instructions  normally  required  by  the  CDC  for  a  job,  the  job 
control  list  (JCL),  the  first  six  lines  of  the  following  code 
The  fourth,  fifth,  and  sixth  lines  reflect  differences  from 
a  general  case  for  program  input  that  are  required  for  the 
specific  case  of  this  simulation  model.  BAMRV  is  the  name 
of  the  permanent  file  containing  the  compiled  FORTRAN  sub¬ 
routines.  The  permanent  file  identification  is  BUZZ.  The 
program,  SLAMPROC,  is  a  procedure  file,  written  in  CDC  con¬ 
trol  language,  that  contains  instructions  for  attaching  and 
processing  the  SLAM  support  program.  The  sixth  line  is  the 
Instruction  for  the  computer  to  begin  the  SLAM  procedure 
file,  using  an  externally  compiled  program,  BAMRV,  as  an 
additional  input,  and  having  a  print  line  limit  of  10,000 
lines. 

The  block  of  code  following  the  JCL  begins  the  SLAM 
portion  of  the  code.  At  the  end  of  the  first  line  of  code 
in  this  block  is  a  number  indicating  the  number  of  runs  the 
user  wishes  to  make.  The  resource  statements,  the  last  two 
lines  of  code  in  this  block,  reflect  the  number  of  endo-inter 
ceptors  and  point  interceptors  that  make  up  the  defense.  The 
remaining  network  code,  in  the  order  discussed  in  the  pre¬ 
vious  five  sections,  would  follow  this  block  of  code.  Th_ 
lines  of  code  requiring  user  change  in  this  block  are  as 
follows 


GEN , WILCRIT , THESIS , 18 / 15/ 19  8 1 , a ; 

RESOURCE /ENI (  b ) ,  1 , 

RESOURCE/ PTI (  c ) ,  2 ; 

In  the  first  line,  the  user  enters  the  number  of  runs  desired 
for  a.  In  the  second  line,  b  indicates  the  desired  number 
of  endo-interceptors .  The  c  in  the  third  line  indicates  the 
number  of  point  interceptors  desired. 

The  last  block  of  code  is  the  initialization  state¬ 
ments  that  set  the  values  of  the  global  variables  in  the 
model.  The  use  of  global  variables  provides  extensive  flex¬ 
ibility  that  permits  the  simulation  of  both  widely  divergent 
RV/BMD  systems  and  minor  variations  on  a  single  system.  As 
a  result  of  the  global  variables,  the  structure  of  the  sys¬ 
tem  network  can  be  altered  without  changing  the  actual  pro¬ 
gramming  of  the  network.  Global  variables  are  used  to  open 
and  close  pathways,  activate  and  deactivate  defenses,  and 
define  the  capabilities  and  vulnerabilities  of  the  system. 
Furthermore,  they  are  used  to  reflect  variations  in  the  sys¬ 
tem  parameters  between  simulation  runs.  Tne  user  must  set 
all  of  these  global  variables  to  specific  values  prior  to 
exercising  the  model.  The  access  to  classified  information 
is  critical  in  setting  appropriate  values  for  these  variables. 


Annotated  FORTRAN  User  Function  Code 


FUNCTION  JSERFUFN) 

CC?1«ON/SCO!1i/ATF!8{l^!.OO(ig«).DDL(!03)iBTNQW.:iiKFA»NSTOPl 
+  NCLNR f  NCRDR > NPRNT » NNRUN » NNSET , NTAPE . 33 ( i 33! » SSL (1 30 ) . 

+  TNEXT « TNOUt  XX ( 1331 

Ce«H0N/RW/X»HNU»Al»Cl»C2»HNUT.AIRDEN»TE»TP 

REAL  HNUU3)  ,  A1  (32)  ,C1  (32) .C2(32)  >HNUT(32)  ,«UTR0»«FPr 
+  AIRBENI87) iINTALEilNTAL? 

C 

C  VARIABLE  KEY 

C  X  :  REENTRY  ANGLE  OF  AN  RV  IN  DEGREES 

C  HYP  =  SLANT  RANGE  THAT  RV  FLIES  THROUGH  ENDO  OR  POINT 

C  DEFENSE  ZONES  IN  FEET 

C  IE  --  TIME  FOR  RV  TO  FLY  THROUGH  ENDO  ZONE  IN  SECONDS 

C  UE  =  MEAN  TIME  OF  ABN  INTERCEPT  IN  ENDO  ZONE  IN  SECONDS 

C  SBE  =  STANDARD  DEVIATION  OF  TIME  OF  ABM  INTERCEPT  IN 
C  ENDO  ZONE  IN  SECONDS 

C  RN  =  NORMAL  RANDOM  VARIATE 

C  TP  =  TINE  FOR  RV  TO  FLY  THROUGH  POINT  ZONE  IN  SECONDS 

C  UF  -  KEAN  TIME  OF  ABM  INTERCEPT  IN  POINT  ZONE  IN  SECONDS 

C  SDP  =  STANDARD  DEVIATION  OF  TIKE  OF  ABN  !NTEKC£PT  IN 
C  POINT  ZONE  IN  SECONDS 

C  YIELD  =  YIELD  OF  AN  ASK  OR  RV  IN  KILOTONS 
C  TEMP  -  TEMPERATURE  OF  X-RAYS  IN  KEV 

C  SR  =  MISS  DISTANCE  OF  AN  ABM  OR  RV  IN  FEET  OR  KILOMETERS 

C  DNSS  =  SURE  SAFE  NEUTRON  INTENSITY  IN  NEUTR0NS/CN2 

C  DNSK  =  SURE  KILL  NEUTRON  INTENSITY  IN  NEUTRONS/CHZ 

C  DXSS  -  SURE  SAFE  X-RAY  INTENSITY  IN  CAL0RIES/CM2 

C  DXSK  *  SURE  KILL  X-RAY  INTENSITY  IN  CAL0RIES/CM2 

C  R  =  MISS  DISTANCE  IN  CENTIMETERS 

C  INTALE  =  ALTITUDE  3F  ENDO  INTERCEPT  IN  FEET  OR  KILOMETERS 

C  INTALP  =  ALTITUDE  GF  POINT  INTERCEPT  IN  FEET  OR  KILOMETERS 
C  l  =  ALTITUDE  OF  ENDO  OR  POINT  INTERCEPT  TT  NEAREST  KILOMETER 

C  PM I  =  MASS  INTEGRAL  IN  GM/CM2  (AIR  DENSITY  TIKES  SLANT  RANGE) 

C  PIR2F  *  NEUTRON  FLUENCE  AT  A  PARTICULAR  RANGE  SR  AND  ALTITUDE  L 
C  F  =  NEUTRON  FLUENCE  ON  RV 

C  YIELD!  =  YIELD  OF  ENDO  INTERCEPTOR  IN  CALORIES 
C  SUMF  -  TOTAL  X-RAY  FREE  FIELD  FLUENCE  ON  RV 
C  HNU1  -  NAN-NORMALIZED  ENERGY  OF  X-RAYS  IN  A  PARTICULAR  GROUP 
C  MUTRO  =  X-RAY  ATTENUATION  COEFFICIENT  FOR  AIR  IN  CH2/GK 

C  MFP  =  MEAN  FREE  PATHS  (MASS  INTEGRAL  TIMES  X-RAY  ATTENUATION 

C  COEFFICIENT  FOR  AIR! 

C  BUF  =  X-RAY  BUILD  UP  FACTOR 
C  FFFG  =  FREE  FIELB  X-RAY  FLUENCE  FOR  A  PARTICULAR  GROUP 

C  FFFT  -  TOTAL  X-RAY  FLUENCE  ON  RV 

C  RNSK  =  RANGE  GF  NEUTRON  SURE  KILL  IN  KILOMETERS  OR  CENTIMETERS 

C  RNSS  -  RANGE  OF  NEUTRON  SURE  SAFE  IN  KILOMETERS  OR  CENTIMETERS 

C  DNA  -  ALPHA  VALUE  FOR  NEUTRON  CALCULATIONS 
C  DNB  =  BETA  VALUE  FOR  NEUTRON  CALCULATIONS 
C  UR  OR  URi  -  WEAPON  RADIUS  FOR  NEUTRONS 
C  PR08A  OR  PROB1A  =  POWER  OF  EXPONENT  IS  SINGLE  SHOT  PROBABILITY 


OF  NEUTRON  KILL  CALCULATIONS 

PROS  OR  PROB1  =  SINGLE  SHOT  PROBABILITT  OF  KILL  FOR  NEUTRONS  ' 

RISK  *  RANGE  OF  X-RAY  SURE  KILL  IN  KILOMETERS  OR  CENTIMETERS 
RISS  -  RANGE  OF  X-RAY  SURE  SAFE  iH  KILOMETERS  OR  CENTIMETERS 
DXA  =  ALPHA  VALUE  FOR  X-RAY  CALCULATIONS 
DIB  =  BETA  VALUE  FOR  X-RAT  CALCULATIONS 
NR2  =  WEAPON  RADIUS  FOR  X-RAYS 
PR082A  -  POWER  OF  EXPONENT  IN  SINGLE  SHOT  PROBABILITT  OF 
X-RAY  KILL  CALCULATIONS 

PR0B2  :  SINGLE  SHOT  PROBABILITT  OF  KILL  FOR  X-RATS 

PS  =  SINGLE  SHOT  PROBABILITT  OF  SURVIVAL  OF  RV  IN  ENDO  ZONE 

DPSS  =  SURE  SAFE  OVERPRESSURE  FOR  TARGET 

DPSK  =  SURE  KILL  OVERPRESSURE  FOR  TARGET 

DPA  =  ALPHA  VALUE  FOR  OVERPRESSURE  CALCULATIONS 

BPB  =  BETA  VALUE  FOR  OVERPRESSURE  CALCULATIONS 

SSR  =  MISS  INSTANCE  OF  RV  SCALED  TO  1  KILOTON  IN  FEET 

RR1  =  SSR 

STMRM  =  SCALED  HISS  DISTANCE  OF  RV  IN  METERS 
0P=  OVERPRESSURE  ON  TARGET 

ANGTTi  P9»  =  USED  IN  EMPIRICAL  FIT  FOR  OVERPRESSURE  CALCULATIONS 

R1  :  SCALED  MISS  DISTANCE  OF  RV  IN  KILOMETERS 

DPI  =  ALTERNATE  OVERPRESSURE  FIGURE  ON  TARGET 

Z  =  PROBABILITT  OF  DAMAGE  TO  TARGET 

11  --  ALTERNATE  PROBABILITT  OF  DAMAGE  FIGURE  TO  TARGET 

FROB  =  AVERAGE  TARGET  DESTRUCTION  TO  TARGET  AREA 


INPUT  DATS 


MIDPOINTS  OF  18  EQUAL  ENERGY  GROUPS  OF  PLANCK'S  SPECTRUM 
DATA  KNU/. 77, 1.82, 2.34, 2. 81, 3.27, 3. 76, 4.32,5. 00, 5. 97, 8. 28/ 

THE  NEXT  FOUR  DATA  GROUPS  ARE  ARRAYS  OF  KALANSKT’S  CONSTANTS 

TO  COMPUTE  SUILD  UP  FACTOR 

DATA  KNUT/ 12. 8, 14.0, 16.0, 18.3,28.0, 

+  22. 3(24. 3(26. 3(23.3(33. 3r 

+  32. 0(34. 3(36.0(33. 0(40. 0t 

+  45. 3(58. 0?S5.0(60. 0(73. 3( 

+  30. 8( 180. 8(123. 8(150. 3(280. 3( 

+  258. 3(388.3(350. 3(488. 3(583. 3( 

*  638.8(750.3/ 

DATA  Al/-. 227, -.373, -.323, -.634, -1.072, 

♦  -1.048,-1.7*8,-2.673,-2.664,-6.338, 

+  -3.805,-8.504,-75.833,-20.330,-16.940, 

+  14. 59(11. 31(109.2(-U.35(-114.1( 

+  -113.1,-10.93,-8.153,13.14,-38.92, 

4  -6,388(19. 89(72.37, -6. 063(16.29( 

+  -57.58,-17.2/ 


DATA  Cl/-.4(-.4,-.68(-.46»-.36( 
+  -.48, -.34, -.26, -.3,-. 14, 

+  -.1,-. 1,8.0,-. 02, -.02, 

+  .12, .16, .12,0.0, .14, 

+  .16, .36, .02, .26, .18, 

+  3.0, .2, .16,0.8, .16, 

+  .1,.06/ 


DATA  C2/8. 0,0.0, .32, .02, .02, 
+  .84, .84, .84, .06, .04, 

+  .04,. 86,. 82,. 36, .08, 

+  -.02»-.04,.l,.2,.16i 

+  . 18, .26, .28, .1, .2, 

4  .26, .12, .14,. 22, .08, 

+  .12,. 12/ 


AIR  DENSITY  IN  KC/H3  FOR  ALTITUDE  FROM  0  TO  86  KM 
DATA  AIRDEN/1. 225, 1.1117, 1.0066, .90925,. 81935,. 73643, .66011, 
+  .59002, .52579,. 46706,. 41351,. 36480,  1194, .26660, 

+  .22786, .19476, .16647, .14230,. 12165  .8400, .88891, 

+  .075715,. 06451,. 055006,. 046938,. 04  ^84, .034257,. 029298, 

+  .025076,. 021478, .018410, .015792, .013555, .011573, 

4  . 0098874, , 0084634, . 0072579 , . 0062355, . 0053666, . 0046268, 

4  .0339957, .0034564, .0029948, .8025989, .8022589, .3019663, 

4  .8017142, .0014965, .0013167, .8011628, .0010269, .0009069, 

4  .00080562,. 80071791,. 00063901 (.00056810,. 00050445, 

4  .00044738, .00039627,. 00035054, .00830968, .00027321, 
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0  DETERNINE  RV  REEHTRT  ANGIE  ST  RV  TYPE 
If  (A'SIB(2>  .E  a.m*XI(£J 
If  lATRIBttJ.GT.UMIttJ 
53  7Q  (Ii2i'5ii)5i6i7i3) » IFM 
C 

C  USER  FUNCTIONS  ONE  THROUGH  FIVE 
C 

C  USER  FUNCTION  1  -  CALCULATE  TIKE  TO  FLY  THROUGH  ENDO  DEFENSE  ZONE 
1  HTP  !  80040/S1NDU) 

TE  -  HYP/XX (7) 

USESF  =  TE 
RETURN 
C 

C  USER  FUNCTION  l  -  CALCULATE  TINE  OF  INTERCEPT  IN  ENOO  DEFENSE  ZONE 
2  UE  -  (TE-.D/2 
30E  =  (TE  -  UE)/3 
RN  -  RNGRM(UE,S0£,21 
IHRN.LT. 4. DUSERF^.l 
IF(RN.GT.TE)USERF=TE 
IFHRN.GEJ.  1)  .AND.  (RN.LE.TE)  )US£RF:RN 
RETURN 


C**m**4tlt***Mm4mm4tm*tt*m*4*mm«*4tKm*m**»4U*t 

■  C  USER  FUNCTION  3  -  CALCULATE  TINE  TO  FLY  THROUGH  POINT  DEFENSE  ZONE 
3  HTP  =  Mflf/SINOM 
TP  -  HTP/XX07) 

USERF  *  TP 
RETURN 
C 

C  USER  FUNCTION  4  -  CALCULATE  TINE  OF  INTERCEPT  IN  POINT  DEFENSE  ZONE 
4  UP  =  (TP  -  .i)/2 
SOP  -  (TP  -  UP)/3 
RN  -  RN0RH(UP»SDP,3! 

IFIRN.LT. 4.1)USERF=.l 
IF (RN.GT ,TP)USERF=TP 
IFURN.GE.i.l)  .AND.  (RN.LE.TPi  )U8ERF=RM 
RETURN 


C  USER  FUNCTION  5  -  CALCULATE  TINE  TO  FLT  LAST  13,444  FEET  TO  TARGET 
5  HTP  :  10141/SIND (I) 

USERF  -  HYP/XX (7) 

RETURN 

C 


saEasisseisassMBaefi^, 
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C  USES  FUNCTION  b  -  CALCULATE  SINGLE  SNOT  KILL  PROBABILITY  OF 
C  ENOO  INTERCEPTOR 

4  HELD  =  XXI17) 

TSHP  *  IX (91 
RN=ATRI8(8) 

C  CALCULATE  KISS  DISTANCE 

5R=SQRT  ( XI  ( 1 8)  «2*  (LOG !  1-RN)  /  (-.643) )  i 
PRINT  ♦»'  RANGE  =  * .SR 
IF (ATRIBI2) .EQ. 1)TH£N 

C  SET  X-RAY  AND  NEUTRON  KILL  LEVELS  FOR  ARV 

DNSS  -  11(12) 

DNSK  =  11(14) 

DXSS  =  IX(27) 

DISK  =  IIC2S) 

ELSE 

C  SET  X-RAT  AND  NEUTRON  KILL  LEVELS  FOR  DRV 

DNSS  *  11(29) 

DNSK  =  IK3D) 

DXSS  =  XX (31) 

DXSK  s  XX (32) 

ENDIF 

C  CALCULATE  KISS  DISTANCE  IN  KH 
SR  =  SR  *  .0003048 

C  CALCULATE  KISS  DISTANCE  IN  CK 

R  *  SR  4  If,.  44  5. 

C  CALCULATE  ALTITUDE  OF  ENDO  INTERCEPT  TO  NEAREST  KK 
INTALE  -  150004'.  -  ATRIB(4)  4  XX(7)  4  SIND(X) 

INTALE  =  INTALE  *  .0003048 
L  =  NINTt INTALE)  +  1 

C  CALCULATE  KASS  INTEGRAL 

PHI  *  AIRDEN(L)  *  SR  4  10  44  2 

C  CALCULATE  NEUTRON  FLUENCE  ON  RV 

PIR2F  =  EXP(-6.775  +  .005249  4  PK!  -  .000005344 
+  *  PHI  it  2  -  .00021448  4  PHI  44  i.5  - 

♦  3.8214  4  PHI  4i  .5  +  10.875  4  PH!  44 

+  (1./3.)  -  1.3975  ♦  LOG (PKI)5 

F  =  (3.14E23  *  YIELD  *  PIR2F!/(4*3,1414*fi«2) 

C 

C  CALCULATE  X-RAT  FLUENCE  ON  RV 

C 

C  CALCULATE  HELD  IN  CALORIES 

YIELD!  =  HELD  *  IEI2 

SUHF  =  0. 

DO  101  1  :  1.10 

C  CALCULATE  ENERGT  OF  X-RATS  IN  CROUP  I 

HNU1  =  TEHP  *  HNU(I) 

C  CALCULATE  X-RAT  ATTENUATION  COEFFICIENT  FOR  AIR 

HUTRO  *  -.001354  +  19.7544  /  HNU1  -  441.743 
+  /  HNU1  «  2.  +  4480.02  /  HNU1  »  3, 


A  -  Ml.lt>  I  mi  «  4.  +  937.358  /  RWJ1  »+  5. 

C  LOCATE  KALANSXT'S  CONSTANTS 
HFP  *  PHI  ♦  MUTRO 

IF  tHN'Jl  .LT.  12. )  THEM 
BUF  --  i 
CO  TO  20 
ELSE 

00  230  J  -  h32 
IFIHNUT(.i)  .GT.HNUDTHEN 

IF ( HNUT ( J) -HNU 1 . GT . HNU 1 - HNUT 'J- 1 1 ) THEN 
K  =  J  -  l 
ELSE 

K  =  J 

ENDIF 
GO  TO  10 
END  IF 

20#  CONTINUE 
ENDIF 

C  CALCUALTE  BUILD  UP  FACTOR 

10  A2  -  1  -  AUK) 

BUF  =  AUK)  *  EXPICKK)  *  RFP)  +  A2  *  EXP  (C2 
+  (K)  *  HFP) 

C  CALCULATE  ENERGY  IN  GROUP  I 

20  FFFG  =  ,1  *  BUF  *  EXP(-HUTRO  *  PR 
SURF  =  SURF  +  FFFG 
100  CONTINUE 

C  CALCULATE  TOTAL  X-RAY  FLUENCE  ON  RV 

FFFT  =  .8  ♦  YIELD1  *  ($URF/(4  *  3.1416  *  R  «  2.1) 

PRINT  »(’  F  -  ’tFt'  PIR2F  =  ' rPIRZF 
PRINT  V  FFFT  =  »»FFFT»»  SURF  :'i SURF 
C 

C  ITERATE  TO  FIND  NEUTRON  SURE  KILL  RANGE 
DO  500  JJ= W103 

PHI  =  AIROEN(L)*JJ+10**2 
PIR2F=EXP (-6. 775+. 005269+PRI-. 000015346 
+  +PHIh2-.0#0Z1468»PHI*H.5- 

+  3,8Z14*PHI»*.5+10.875*PHI« 

+  (1 . /3. I -1 .3975+LOG (PHI ) ) 

Fs3.16E23+YIELD*PIR2F/{4*3.1416*yj*l0«5)«2) 
IFIF.L£.DNSK)TH£N 
RNSK  -  JJ 

PRINT  F  =  »(F»»  PIR2F  =  »tPIR2F»»  RNSK  =  ’(RNSK 
GO  TO  90# 

ENDIF 
5fi  CONTINUE 
91#  RNSK-RNSK+10++5 
C 

C  ITERATE  TO  FIND  NEUTRON  SURE  SAFE  RANGE 
DO  501  KK:l(100 

PR  I s  A IRDEN  ( L)  *KK+ 1 0«2 
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PIRZF^ESP  (-6.775+ .  305Z&9*PSI-  .300805346 
*  »PMI«2-.§002l4S3iPHI«1.5- 

+  (l./3.M.3975*LGG(PHM 

F  =  3,16E23»nELD*PIR2F/(4*3.14l6*(KK*10**5)H2) 

IF (F.LE. DNSS) THEN 
SNS3-KK 

PRINT  F  =  *tFj*  PIR2F  -  -»PIR2F»'  RHSS  -  * iRNSS 
GO  TO  m 
END  IF 

501  CONTINUE 

901  RHSS  =  RNSS  *  10  »S 

C 

C  CALCULATE  PROBABILITY  OF  NEUTRON  KILL 
DNA  =  .5  *  LOG  (RNSK  *  RNSS) 

ONB  =  (,5/2.054)*LQG(RNSK/RNSS) 

C  CALCULATE  1IEAPGN  RADIUS  FOR  NEUTRONS 
«R1  =  EXP (DNA  +  DNB  **  2.) 

PRINT  *»»  DNA  =  ’iDNAt’  DNB  =  !?DNBt*  HR1  =  ’»HR1 
C  CALCULATE  SINGLE  SHOT  PROBABILITY  OF  NUETRQN  KILL 
PRGB1A  =  -.693+1KR1/ (XX(18)*30.4S))**2. 
IF(ft'.S<?RGBlA).GT.23)THEN 
PROB1  =  1. 

ELSE 

PROB1  =  1-EXP (PROB1A) 

ENDIF 

C 

C  ITERATE  TO  FIND  X-RAT  SURE  KILL  RANGE 
DO  502  LL=lil00 

RR  =  FLOAT (LL5*  ,2 
PNI  =  AIRBEN(L)+RR*10»2. 

SUKF  =  0. 

DO  503  LLI  *  1)10 

HNU1  =  TEHP  *  HNU(LLI) 

HUTRO  =  -. 001354+19. 7564/HNU1-461. 763 
+  /KNUl«2.+6680.02/HNUi«3. 

+  -3497.36/KNUl«4.+907.358/HNUl«5. 

HFP  =  PHI  »  HUTRO 
IFlHNUl.LT. 12. JTHEN 
BUF  =  1. 

GO  TO  902 
ELSE 

DO  504  HR  =  1 »32 

IF(HNUT(HH).GT.HNU1)THEN 

IF(HNUT(HH)-HNU1.GT.HNU1-HNUT(HH-1) JTHEN 
K  =  HH  -  1 
ELSE 

K  *  HH 
ENDIF 
GO  TO  903 
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903  42  -  l-AKK) 

BUF  =  Al(K)*EXP(Cl(K)*NFP)+A2*EXP{C2 
+  (KNHFP) 

902  FFFC  =  .1»8UF«EXP{-HUTR0»PNI) 

SURF  -  SURF  +  FFFC 
503  CONTINUE 

FFFT  =  .8*nELBHSUNF/(4»3.1416*{RR*10M5.)«2.) 

IF (FFFT.LE.DXSKS THEN 
RXSK  =  RR 

PRINT  +,'  FFFT  -  '  .FFFT.'  SURF  -  ».SUNF 
PRINT  RXSK  =  '.RXSK 

GO  TO  904 
ENOIF 
502  CONTINUE 
984  RXSK=RXSK*13+*5 
C 

C  ITERATE  TO  FIND  X-RAT  SURE  SAFE  RANGE 
DO  505  11=1.100 
RR=FLOAT  (LL)  ■» .  2 
PNI  =  flIRDEN(L)*RR*18*+2. 

SURF  "  0 1 
DO  506  ILI  =  1.10 

HNU1  =  TENP  *  ilNU(LLI) 

NUTRO  =  231354+19. 7564/HNUI-461. 763 

+  /HNUl«2.+6680.02/HNUl«3. 

+  -3497 . 36/HNUl«4 .  +907 .358/HNUU+5 . 

NFP  =  PNI  »  NUTRO 
IFIHNU1.LT.12.1THEN 
BUF  =  1. 

CO  TO  905 
ELSE 

DO  507  H«  =  1.32 

IFIHNUTINNLGT.HNUDTHEN 

IF (HNUT INN) -HNU1 .CT .HNUi-HNUT (NN- 1 ) ) THEN 
K  =  NN  -  1 
ELSE 

K  =  NN 
ENDIr 
CO  TO  906 
ENOIF 


507 

CONTINUE 

ENOIF 

906 

A2  =  l-AKK) 

BUF  =  A1 (K) +EXP (Cl (K) +NFPS+A2+EXP (C2 

+ 

(K)*NFP) 

905 

FFFC  =  .UEUFJEXP(-NUTRO*PNI) 

SURF  =  SURF  +  FfFG 

905 


39S 


CONTINUE 

P~T  ;  ,a*YiELDl*SU»Fm*3,14U*(RR*10«5.)«2.) 

IF (f?FT .Lt.DXSS) THEM 
RXSS  -  SR 

PRINT  *,'  FFFT  -  *»FFFTr»  SUflF  =  >t8UHF 
PRINT  *»’  RXSS  *  ',RX$S 
GO  TO  907 
ENDIF 
505  CONTINUE 
907  RXSS=RXSS* 10**5 
C 

C  CALCULATE  PROBABILITY  OF  X-RAT  KILL 
OXA  *  .5*L0G(RX$K*RXS$) 

DXB  =  (.5/2.054) *LOG (RXSX/RXSS! 

C  CALCULATE  WEAPON  RADIUS  FOR  X-RAYS 
WR2  -  EXP(DXA*0X8**2.) 

PRINT*,'  OXA  *  ’iDXA,’  DXB  s  ’,0X8,'  UR2  *  MIR2 
C  CALCULATE  SINGLE  SHOT  PROBABILITY  OF  X-RAT  KILL 
PR0B2A  =  -.693*(WR2/<XX(18)*30.48))**2. 

IF(ABS(PR0B2A) .GT.23.1THEN 
PR0B2  =  1. 

ELSE 

PR0B2  -  1  -  £XP(PRflB2A) 

ENDIF 

C  CALCULATE  SINGLE  SHOT  PROBABILITY  OF  RV  SURVIVAL 
PS  =  (1-PROBD*  (1-PR0B2) 

C  CALCULATE  SINGLE  SHOT  PROBABILITY  OF  RV  KILL  BY  ENDO  INTERCEPTOR 
USERF  =  1  -  PS 


LSER  FUNCTION  7  -  CALCULATE  SINGLE  SHOT  KILL  PROBABILITY  OF 
POINT  INTERCEPTOR 

7  TUILD  -  XXU9) 
fM-ATRIR(ll) 

CALCULATE  MISS  DISTANCE 
SF'-SQKTin  (20)  4*2*  (LOG  ( 1-RN)  /  <- .49311) 

Pfl.'IT  *,'  RANGE  =  ’.SR 
IF.ATRI8(2).Ei.l!THEN 

SET  NEUTRON  KILL  LEVELS  rOR  ARV 
DNSS  -  ■1(121 
BNSK  =  aX (14) 

ELSE 

SET  NEUTRON  KILL  LEVELS  FOR  8RV 
DNSS  =  71(29) 

DNSK  *  >.2(301 
ENOIF 

CALCULATE  MISS  DISTANCE  IN  KM 
SR  s  SR  *  .000304:; 

CALCULATE  MISS  DISTANCE  IN  CM 
R  *  SR  *  10  *>  5 

CALCULATE  ALTITUDE  OF  ''01  NT  INTERCEPT  TO  NEAREST  KM 
INTALP  =  70F00.  -  ATi)IE:4)  *  XX(7)  *  SIND(X) 

INTALP  =  5KTALP  *  .0003048 

M  8  N1NT (INTALP)  +  t 

CALCULATE  MACS  INTEGRAL 

PHI  =  AIRDEN(M)  *  SR  *  10  **  2 

CALCULATE  NE1/TR0N  FLUENCE  ON  RV 

PIR2F  =  EXP (-4. 775  +  .005249  *  PMI  -.000005344 


*  *  PMI  «  2  -.00021463  *  FMI  **  i.S  - 

*  3.8214  <•  PMI  «  .5  +  10,875  *  PMI  H 

*  (1./3.)  -  1.3975  *  LOG  (PHD) 

F  8  (3. 16£23  *  HELD  *  PIR2F17  (4  *  3.1414 
+  *  R  **  2) 

PRINT*,’  F  8  »,fV  PIR2F  8  >,PIR2F 


ITERATE  TO  FIND  NEUTRON  SURE  KILL  RANGE 
DO  508  JJ=t,100 
RR=FL0AT(JJ)*.2 
PMI  8  AIRD£N(M)*RR*13**2 
P IR2F=EXP ( -4. 775* . 005269*PMI - .000005346 
+  *PMI**2-.00021468*PMI«t.5- 

+  3.8214*PMI**.5+10.875*PMI« 

*  (l./3.)-!.3975*L0C(PHI)> 

F=3.16E23*iIELD*PIR2F/(4*3.1414*(RR*10h5)**2) 
IF(F.LE.DNSK)THEN 
RNSK  8  RR 

PRINT*,’  F  =  » ,F»  *  PIR2F  ?  ’,PIR2F,’  RNSK  ==  »tRNSK 
CO  TA  908 
END  IF 


588  CONTINUE 
838  RNSK-RNSK* 18**5 
C 

C  ITERATE  TO  FIND  NEUTRON  SURE  SAFE  RANGE 
CO  539  KKMtill 
RR-FL0AT(KK)*.2 
PHi-AIRDEN(N'*RP.*l.3«2 
PIR2F=£XP(-4.775*.885269*PRI-. 038885344 
♦  *PNI«2-.03321448*PHI**I.5- 

+  3.8214*PRI«.5+10.875*PNI« 

+  U./3.M.3975*L0C(PHIi) 

F  =  3.14E23*nELD*PlR2F/(4*3.1414*(RR*13**5)«2) 

I F ( F . LE . DNSS ) THEN 
RNSS=RR 

PRINT  ♦.»  F  =  »tF»'  PIR2F  =  #»PIR2F.»  RNSS  =  ».RNSS 
CO  TO  969 
END  IF 

589  CONTINUE 

989  RNSS  =  RNSS  *  16  **5 

C 

C  CALCULATE  PROBABILITY  OF  NEUTRON  KILL 
DNA  :  .5  *  LOG  (RNSK  *  RNSS) 

0N8  =  ( .5/2.854) *LGC (RNSK/RNSS) 

C  CALCULATE  WEAPON  RADIUS  FOR  NEUTRONS 
WR  =  EXP(DNA*DNB«2) 

PRINT  DNA  =  1  iDNAi’  SNB  =  ’iDN3»»  WR  *  MIR 
PROBA  *  -,493*{WR/<XX<28)*33.48)!«2. 

C  CALCULATE  SINGLE  SHOT  PROBABILITY  OF  RV  KILL  BT  POINT  INTERCEPTOR 
IF(ABS{PR08A).GT.23.)THEN 
PROB  =  1. 

ELSE 

PROB  =  1  -  EIP(PRQBA) 

ENDIF 

USERF  =  PROB 
RETURN 


04444*44444444444444444*4444*44**4*444444*444444444444444*44444444444 

C  USER  FUNCTION  3  -  PROBABILITY  THAT  SV  DESTROYS  TARGET 
8  RN-ATR!Jil4) 

iF{ft7RIB(2).E8.i)THEN 
HELD  --  11(21) 

C  CALCULATE  HISS  DISTANCE  OF  ARV 

SR=S5RT  IXX  (22)  «2i(L0G(l-RS)  /  (-  .4,93) ) ) 

ELSE 

HELD  =  IMS) 

C  CALCULATE  HISS  DISTANCE  OF  DRV 

SR=SGRT(XX (24)442* (LOG ( 1-RS) / (- .893) ) ) 

END  IF 

PRINT  *»'  RANGE  =  '.SR 
C  SET  SURE  SAFE  OVERPRESSURE  FOR  TARGET 
DPSS  =  IX (S3) 

C  SET  SURE  KILL  OVERPRESSURE  FOR  TARGET 
OPSK  =  XI (34) 

DPA  =  .5  *  LOG  (DPSK  *  DPSS) 

BPS  =  (.5/2.854)  *  LOG  (DPSK/BPSS) 

C  SCALE  HISS  DISTANCE  TO  1  KILOTGN  EQUIVALENT 
SSR  -  SR  /  HELD  4*  U./3.) 

RRi  -  SSR 

C  CALCULATE  HISS  DISTANCE  IN  H 
STHRH  -  RRI  /  3.2815 

C  EMPIRICAL  FIT  FOR  OVERPRESSURE  CALCULATION 
IF(RR1,E8.0)THEN 

ANGTf  =  3.1416  /  2 
ELSE 

ANCTT  =  0. 

ENDIF 

P9  *  .81  4  EXP (40.3  *  STHRH  **  (-.295)) 

P8  =  .881  *  EXP (31 .3  *  STHRH  «  (-.2136)) 

DP  =  P9  -  (P9-P8)  *  KCOS(ANGTT)  )**2! 

PRINT*,’  DP  =  »fDP 

C  ALTERNATE  OVERPRESSURE  CALCULATION  (NOT  USED  IN  HODEL) 

Rl- (SR/ ( (1 .5*1 IELD) ** (1./3.) ) )*. J883348 
DP1-EIP ( . 19* (LOG (Rl ) ) **2- 1 .5* (LOG (Rl ) ) - . I ) 

PRINT  *»’  DPI  -  ’.DPI 

C  CALCULATE  PROBABILITY  OF  DAHAGE  TO  TARGET 
Z  =  (LOG (DP) -DPA) /DPB 
IF (Z.LT.i.)THEN 
Z  -  ABS  (  Z) 

PROB  »  ,5  *  (  1  +  .1V&854  *  Z  +  .  115194  *  Z  «  2 

*  4  .888344  *  Z  **  3  +  .019527  *  Z  **  4)  «  (-4) 

asE 

PROB  •  1  -  .5  *  (1  +  .196854  *  Z  4  .115194  *  Z  *»  2 

+  4  .188344  *  Z  **  3  4  .019527  *  Z  4*  4)  «  (-4) 

ENDIF  , 

C  CALCULATE  AVERAGE  TARGET  DESTRUCTION  FOR  TARGET  AREA 
PROB  =  PROB  /  11(26) 


ID  4 


Conclusions 


The  entire  code  for  the  model  is  entered  in  the  com¬ 
puter  in  the  following  sequence: 

1)  JCL  cards,  first  block  of  code  in  SLAM  support 
code  section 

2)  *EOR 

3)  Initial  SLAM  code,  second  block  of  code  in  SLAM 
support  code  section 

4)  Arriving  RV  segment  of  code 

5)  Endo  defense  segment  of  code 

6)  Point  defense  swicch  segment  of  code 

7)  Point  defense  segment  of  code 

8)  Target  area  segment  of  code 

9)  Initialize  cards,  last  block  of  code  in  SLAM 
support  code  section 

10)  SLAM  supplementary  code  as  discussed  in  Chapter 
IV  to  carry  out  the  user’s  experiment 

In  order  to  exercise  the  model  from  an  interactive  terminal, 
the  FORTRAN  subroutines  are  then  entered,  compiled  with  a 
file  name,  and  the  resulting  binary  file  is  stored  in  a 
permanent  file  under  the  compiled  file  name.  The  recommended 
method  of  exercising  the  model  is  with  the  interactive  termi¬ 
nal,  as  described.  Card  decks  may  also  be  used.  Using  this 
procedure,  the  FORTRAN  subroutines  would  be  entered  between 
steps  2  and  3  of  the  sequence  above,  followed  by  a  *E0R 
card.  The  user  is  cautioned  to  reference  the  JCL 
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instructions  for  the  particular  computer  when  developing  the 
JCL  cards  or  code. 


APPENDIX  D 

PROGRAM  VERIFICATION 


Introduction 

The  purpose  of  this  appendix  Is  to  illustrate  the 
verification  process  that  ensured  the  computer  model  operated 
as  intended.  The  example  verification  compares  the  output 
of  the  supporting  subroutines  with  the  results  of  analytical 
solutions.  These  analytical  techniques  were  developed  in  the 
three  course  sequence  dealing  with  the  effects  of  nuclear 
weapons  in  the  Strategic  and  Tactical  Science  program  and  are 
based  primarily  on  the  information  contained  in  The  Effects 
of  Nuclear  Weapons  by  Glasstone  and  Dolan  and  A  Short  Course 
in  Nuclear  Weapon  Effects  by  Bridgman  and  John. 

The  verification  shows  a  comparison  of  results  of  an 
engagement  in  the  endo-interceptor  region  using  neutron  flu- 
ence.  This  is  not  a  test  of  the  random  variate  generation 
routines  of  either  the  IMSL  or  the  SLAM  program.  Those  are 
assumed  as  valid  for  the  purposes  of  this  appendix. 

Engagement  Verification 

The  attacker  is  a  maneuvering  RV  with  a  yield  of  150 
kilotons  and  a  CEP  of  200  feet.  The  hardness  levels  of  the 
RV  are  shown  in  Table  XV.  The  defending  ABM  has  a  yield  of 
1  megaton  and  a  CEP  of  3000  feet. 

The  variables,  DNA  and  DNB,  can  be  defined  by  Eqs 
(56)  and  (58). 

DNA  =  .5*  ln(RNSK  *  RNSS)  (105) 

MB  -  -27 hr  * ln  -Hi-  <106> 
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TABLE  XV 


MRV  Hardness  Levels 

X-Rays (cal/cm2 )  Neutrons (neutrons /cm2 ) 


Sure 

Safe 

70 

1  X 

10 13 

Sure 

Kill 

100 

1  X 

1013 

RNSK  and  RNSS  are  dependent  upon  the  radiation  levels  of  the 
engagement.  Thus,  these  must  be  calculated.  The  neutron 
fluence  calculation  will  be  used  to  illustrate  the  process. 

Assuming  a  miss  distance  of  853.8  feet  and  an  inter¬ 
ception  altitude  of  22.6744  km  (74,391  ft),  then  the  neutron 
calculation  proceeds  as  follows. 

The  mass  integral,  MI,  is 

MI  =  (air  density)  *  slant  range  *  10*  (107) 

MI  *  (.055006)  *  (853.8  *  .0003048)  *  102 

*  1.4315  gm/cm2  (108) 

The  ( 4ttR2 )  fluence  is  determined  using  Eq  (41). 

PIR2F  *  e(“6.775  +  .005269(MI)  -  . 000005346(MI) 2 
-.00021468(MI)1,5  -  3.8214(MI)‘5 
+  10. 875(MI) 1/3  -  1.3975  ln(MI))  (109) 

Since  MI  =  1.4315s  PIR2F  =  1.5144,  the  actual  neutron  fluence 
is,  then. 


150 


P  =  (3.16E23  *  1000  *  1. 51410/(4  *tt*( 853 . 8  *  30. 48)2) 

=  5.6231  x  lO'*  neutrons/cma  (110) 

which  agrees  with  the  output  of  the  program,  5.6231  x  lO1^. 

Since  the  ranges  for  sure  safe  and  sure  kill  fluences, 
RNSS  and  RNSK,  respectively,  are  dependent  of  the  mass  inte¬ 
gral,  they  cannot  be  calculated  directly.  An  iterative 
process  is  used  which  yields  RNSK  =  4km  and  RNSS  *  16km. 
Therefore,  using  Eqs  (56)  and  (58),  DNA  =  13*592  and  DNB  * 
-.337.  Eq  (80)  is  used  to  calculate  the  weapon  radius  (WR). 

HR  =  e(DNA  +  tDNB)2>  =  896,1.66.  (Ill) 


Eq  (103)  results  in  the  probability  of  kill. 
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